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ABSTRACT 


This thesis develops a robust analytical computer simulation model of the 
Extended Close Battlefield (ECB) to examine the performance of the command, control 
and communications (C°) systems of the Extended Fiber Optic Guided Missile (EFOGM). 
using both digitized and non-digitized communications. The ECB is represented by a 24 
state Semi- Markov chain formulation. It contains transient and absorbing states that the 
simulation models using Monte Carlo processes and probablistic time distributions. The 
primary measure of effectiveness (MOE) used for comparison between the digitized and 
non-digitized systems is the total time required to process a call for fire (CFF). Sensitivity 
analysis is performed to compare the amount of time a CFF spends in a queue. waiting to 
be processed, within the non-digitized system over a vanable range of available targets. 
Additional sensitivity analysis 1s performed by adjusting the input time parameters of the 
probablistic time distributions to replicate the stress of continuous combat operations. 
While the amount of time a CFF spends in a queue can be brought to zero for specific 
numbers of available targets, the digitized system outperforms the non-digitized system 
over all nee of available targets. using both standard and increased values of input time 


parameters. 





THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this thesis may not 
have been exercised for all cases of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of computational errors, they cannot 
be considered validated. Any application of these programs without additional verification 
is at the risk of the user. 

Additionally, a portion of the analysis conducted for this thesis was performed 
using APL2/PC and AGSS. Naval Postgraduate School uses this program under a test 


agreement with IBM Research. 
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EXECUTIVE SUMMARY 


The Rapid Force Projection Initiative (RFPI) is a study sponsored by the 
Dismounted Battlespace Battle Laboratory (DBBL) located at Fort Benning, Georgia. 
The capabilities to be demonstrated through the RFPI are those that will enable the 
maneuver Commander to engage enemy forces before the enemy can engage him in the 
Close Battle. Specifically, RFPI is to provide new capabilities within the 5 to 15 kilometer 
band of the battlefield known as the Extended Close Battlefield (ECB). The significance 
of these capabilities to the dismounted soldier is in the improvement of supporting sensors, 
command, control and communications (e) devices, and non-line-of-sight (NLOS) 
weapons. These new capabilities will enable fnendly commanders to significantly reduce 
enemy combat power before U.S. soldiers are involved in the direct fire battle. One of the 
NLOS weapon systems proposed under the RFPI is the Enhanced Fiber Optic Guided 
Missile (EFOGM). 

EFOGM is a fiber optic guided missile capable of engaging and destroying enemy 
weapon systems at extended ranges. The system can engage targets within a full 360 
degree radius from its battery location. The missile is a top attack round that allows the 
gunner to adjust the missile’s trajectory up to the point of impact, and is guided through a 
fiber optic data link. The data link connects the missile’s seeker head to a console at the 
gunner’s station. This NLOS weapon greatly improves the ability of ground force 
commanders to influence the ECB by destroying high value enemy targets before the 
enemy has a chance to execute a coordinated, synchronized attack. 

Since EFOGM is a NLOS system, it must depend primarily upon forward placed 
observers and sensors to provide target information. The EFOGM targeting system will 
be integrated into the Advanced Field Artillery Tactical Data System (AFATDS), which 
will allow rapid transmission of target data from the observer/sensor to the battery using 
digitized communications. This will permit parallel processing of targeting information 


versus the current series processing inherent in voice communications. It is logical to 
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assume that the expeditious delivery of target information would enhance the ability of the 
EFOGM to prosecute targets. With a digitized C° system, a sensor or scout could 
simultaneously burst transmit target information, using pre-formatted reports, to all 
headquarters requiring target information, as well as the firing battery. The data could 
then be simultaneously processed for firing at all locations. Current voice C° systems 
require the observer/sensor report to be relayed through a series of single processing 
nodes. Each node must quickly process target information and pass it to the next node in 
the series, until it ultimately becomes firing data as an approved mission at the battery 
level. 

The objective of this thesis is to determine the impact of a digitized C’ system 
versus a non-digitized C’ system concerning target prosecution with the EFOGM. 
Specifically, it addresses the question of whether a digitized communication system 1s 
more effective in handling EFOGM targeting and target prosecution than a voice system. 
This study uses a computer simulation model to examine the performance characteristics 
of the voice and digitized fire control network for the EFOGM. The approach used in 
building the computer model ts to take current procedures for calling indirect fires on 
targets, and apply them to the employment of the EFOGM within the context of the 
computer simulation model. 

By executing a large number of replications of the simulation model, values for 
EFOGM call-for-fire (CFF) processing times and the time those CFFs waited for 
processing are collected. With data from simulation runs of both the digitized and non- 
digitized C’ systems, operating ranges for the two systems are obtained. For each run, the 
specific time events occur are recorded for use in data analysis. These events are included 
for both the digitized and non-digitized replications, and each is recorded for specific 
targets and observers. To compare and contrast the differences between the digitized and 
non-digitized systems, total CFF processing times are calculated by taking the difference 
between the time a target was detected and the time the EFOGM battery finishes 


processing the CFF. Additionally, to test the hypothesis that CFF processing times would 
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increase under the stress of continuous combat operations, values for input time 
parameters are increased by 25% and the replications executed again. 

Examining the data from the output of the simulation replications, 1t was 
discovered that the digitized system is superior to the non-digitized regardless of the 
number of targets presented in a replication. Additionally, if the number of missiles 
available to the battery 1s increased, re-shoots of targets are more prevalent. These 
additional missiles cause an increase in the number of CFFs processed by both the 
digitized and non-digitized C’ systems. While the increase in numbers of CFFs processed 
has a significant effect on the CFF processing times for the non-digitized system, only a 
slight increase was noticed in the CFF processing times for the digitized system. Finally, 
the stress of continuous combat operations was modeled through a 25% increase in 
selected processing times. While this caused a disproportionate increase in the non- 
digitized CFF processing times. processing times for the digitized system remained 
relatively proportional to the 25% increase in selected times. 

This thesis demonstrates the potential of relatively simple, computer simulation 
models to estimate performance parameters in undeveloped combat systems. This 
information can provide analytic agencies with the ability to incorporate future 
technologies into existing combat models, or be used alone for data analysis. Ultimately, 
it provides valuable information for defense planers to preparation for future operations 


and missions around the world. 
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I. INTRODUCTION 


A. BACKGROUND 


Changes in the Department of Defense (DOD) caused by the lessening of tensions 
in Eastern Europe and the lessons learned from Operations Just Cause and Desert Storm 
(ODS) have resulted in DOD budgets being significantly smaller than previous years. 
Additionally, the lessons of ODS have shown that the early deployable land forces need 
more capable and lethal equipment, and lighter and more easily deployable medium to 
heavy force equipment. As a result, the DOD’s Science and Technology (S&T) program 
was established to affordably produce and field new military systems that are more robust, 
but also provide more “bang for the buck”’. To support these goals and to prioritize the 
S&T program, the’DOD established seven S&T “Thrusts” so the highest priority current 
deficiencies could be fixed. Specifically, Thrust Five, Advanced Land Combat, is focused 
to provide new technology for land forces in close combat. Thrust Five consists of three 


top level demonstrations of technology and hardware specific to: 


e Advanced Vehicle Technologies 
e Rapid Force Projection Initiative (RFPI) 
e 21st Century Land Warrior 


Specifically, the goal of the RFPI is to investigate technologies concerned with lethality, 
survivability and deployability to provide improved warfighting capabilities for a force 


projection Army. 
B. RFPI 


The Rapid Force Projection Initiative (RFPI) is sponsored by the Dismounted 
Battlespace Battle Laboratory (DBBL), Fort Benning, Georgia. The capabilities to be 
demonstrated through the RFPI are those that will enable the maneuver commander to 


engage enemy forces before the enemy can engage him in the Close Battlefield. (The 


Close Battlefield 1s that portion of the battlefield reaching from zero to five kilometers in 
front of the foxhole.) Specifically, these initiatives are to provide new capabilities within 
the 5 to 15 kilometer band of the battlefield known as the Extended Close Battlefield 
(ECB). This is especially critical when considering the composition of our light, highly 
deployable forces. The significance of RFPI to the dismounted soldier is in the 
improvement of supporting sensors, command control and communications (C’) devices, 
and non-line-of-sight (NLOS) weapons that will enable friendly commanders to 
significantly reduce enemy combat power before his soldiers are involved in the direct fire 
battle. One of the NLOS weapon systems proposed under the RFPI is the Enhanced Fiber 
Optic Guided Missile (EFOGM). 


C: BASIC DESCRIPTION OF EFOGM 


EFOGM is a fiber optic guided missile capable of engaging and destroying 
camouflaged, dug-in, moving or defiladed targets at extended ranges. The system can 
engage targets within a full 360 degree radius from its battery location. The missile is a top 
attack round that allows the gunner to adjust the missile’s trajectory up to the point of 
impact, and is guided through a fiber optic data link. The data link connects the missile’s 
seeker head to a console at the gunner’s station, which also allows the gunner to scan the 
target area and select targets. This NLOS weapon greatly improves the ability of ground 
force commanders to influence the ECB by destroying high value enemy targets before the 


enemy has a chance to execute a coordinated, cohesive attack. 
1. Employment 


Since EFOGM is a NLOS system, it must depend primarily upon forward placed 
observers/sensors to provide target information. The EFOGM targeting system will be 
integrated into the Advanced Field Artillery Tactical Data System (AFATDS) allowing for 
rapid transmission of target data from observer/sensor to shooter through the use of 
digitized communications. The target data are rapidly processed and sent to the gunner’s 


console, where the gunner will select a route for the missile to fly to the target. Once 


firing confirmation is received, the gunner will launch and track the flight plan until the 
missile reaches the target area. At that point, the gunner can manually guide the missile on 
to the target, or “lock-on” the target and take control of another missile. Some of the 


characteristics of the EFOGM are given in Table 1. 


Maximum Missile Range 15 Kulometers 


Minimum Missile Range 1 Kilometer 


Time of Flight 10 seconds per Kilometer (100 m/s) 
2 missiles within 30 seconds 
Number of On-board Missiles 


Mobility Equal to a High Mobility Multi-Purpose 
Wheeled Vehicle (HMMWV) 


Sustainabilit 96 hours continuous 


Two men 
Table 1. EFOGM Required Characteristics 
jap Targets 


Targets for the EFOGM can be both pre-planned and targets of opportunity. Due 
to the long range and precision of the EFOGM, it will be possible to engage a variety of 
high value targets while the enemy 1s still in march formation, or before moving into the 
effective range of his weapons during an attack. Examples of high value targets are enemy 
command and control assets (C’), air defense assets (ADA), engineering assets (ENG), 
artillery pieces and prime movers, helicopters and selected armored vehicles. Since there 
are only a limited number of EFOGMs available within a rapidly deploying force, it will 


NOT be used as a mass armor killing system. 
3S: Command and Control (C*) 


The EFOGM 1s a divisional asset that is fought at the brigade level. The EFOGM 
company will consist of three platoons of four squads. Target acquisition is provided 
through the Fire Support Element (FSE) located with the brigade tactical operations 
center (TOC) to an EFOGM liaison officer (LNO). The EFOGM liaison element consists 


of a section leader (normally the EFOGM battery commander) and two target analysts. 


The fires of the EFOGM company are coordinated by the brigade Fire Support Officer 
(FSO), who passes targeting information and engagement orders through the EFOGM 
liaison element. The EFOGM LNO then passes firing orders to the EFOGM battery. 
EFOGM platoon headquarters (HQ) directs which squad will execute the fire mission to 
reduce the chance of target duplication. With the introduction of digitized 
communications, targeting information can be sent simultaneously to the brigade FSO, 


EFOGM LNO and the firing battery to reduce engagement time. 
4. Target Acquisition 


Assets that acquire targets for the EFOGM range from those at the national level 
to those in front line combat units. These could include human intelligence (HUMINT) 
reports generated from units in contact, reports from such platforms as the Unmanned 
Aerial Vehicles (UAV), the future scout vehicle with enhanced sensor arrays, or ground 
surveillance radars (GSR). Additionally, the EFOGM can be used to gather intelligence 


and information using its video seeker head while enroute to targets. 


Si. Target Processing and Allocation 


All target information, regardless of source, is ultimately processed at the brigade 
TOC by the brigade FSO. Allocation of EFOGMSs to targets is made in support of the 
brigade commander’s intent and is coordinated through the EFOGM LNO in the TOC, 
who prepares and coordinates EFOGM battery fires. 


D. DIGITIZED COMMUNICATIONS 


There has been much attention focused on the potential benefits of digital 
communication systems on the battlefield. Digitized communication permits parallel 
processing of targeting information versus the current series processing inherent in voice 
communications. It is logical to assume that the expeditious delivery of target information 


would enhance the ability of the EFOGM to prosecute targets. Likewise, any delays in 


transmitting target information could increase the possibility that the target will not be in 
the field-of-view of the seeker when the missile reaches the target area. 

With a digitized communication system, a sensor or Scout can burst transmit target 
information through a pre-formatted Size-Activity-Location-Uniform-Time-Equipment 
(SALUTE) report. The assistance of Global Positioning Systems (GPS), laser range 
finders, and improved frequency hopping Single Channel Ground and Airborne Radios 
(SINCGARS) enhances the accuracy and delivery potential of these forward located 
observer/sensors. As mentioned, information can arrive simultaneously at the division 
ready brigade (DRB) TOC, EFOGM LNO, and the firing battery location. The data 
would automatically be processed for firing at all locations. Using a “do not fire by 
exception” policy, the target will be fired unless told not to by the DRB HQ. Current 
voice C° systems require the observer/sensor report to be relayed through a series of single 
nodes, which must quickly process the information and pass it on to the next node in the 
series, before it ultimately becomes firing data in an approved mission at the firing unit 
level. Handling of multiple or simultaneous fire missions is not possible, so fire missions 
may have to wait valuable time to be processed at each level of control, thereby reducing 


the number of targets that may be engaged in a time period. 
| De STATEMENT OF THE PROBLEM 


There is little information available concerning the benefits of digitized 
communication in terms of normal measures of combat effectiveness. Current high 
resolution wargaming models such as JANUS and CASTFOREM do not adequately 
model delays caused by the processing of calls for fire (CFF) generated by forward 
observers or scouts. The objective of this thesis is to determine the impact of a digitized 
versus non-digitized communication system with respect to target prosecution with the 
EFOGM. Specifically, how much more effective is a digitized communication system in 
handling EFOGM targeting and target prosecution versus the currently available voice 


system? These results may help determine whether a digitized communication system, 


coupled with the EFOGM, will significantly improve the lethality and survivability of 
deployable forces; specifically, the division ready brigade (DRB). Intuitively, one would 
expect that a decrease in the target prosecution time would result in more missiles fired 
over a similar period of time. Additionally, the chances of a EFOGM being "successful" 
(i.e., striking a desired target) should improve because of the increased time the EFOGM 


has to acquire a target in the engagement area. 


F. SCOPE 


The DBBL has developed a high resolution scenario (HRS) in JANUS that 
replicates a DRB in a forced entry operation in Latin America (HRS 33.5). This HRS is 
the foundation for the friendly and enemy force compositions and scenario used in this 
study. Using a state space representation of the HRS, a simple, highly robust computer 
model has been developed to provide a range of values for target prosecution times. 
These values can be input into current high resolution simulation models to replicate the 
delays caused by call for fire (CFF) processing. 

This study will use a discrete time analytical combat simulation model to examine 
the performance characteristics of the voice and digitized fire control network for the 
EFOGM missile. The combat environment in which the EFOGM operates will be 
modeled as a Semi-Markov process with transient and absorbing states, combining a 
Monte Carlo process within the simulation to model the events of acquiring, processing, 
and prosecuting targets by the EFOGM. The approach used in building the state space 
description was to take current procedures for calling indirect fire on targets, apply them 
to the employment of the EFOGM and model these steps as states in a state space 


diagram. The description of this state space is shown in the following chapter. 
G. STUDY LIMITATIONS AND ASSUMPTIONS 


While this study will seek to closely replicate engagements by the EFGOM on 


targets in the ECB, there are certain aspects of the ECB that were not, or could not be 


covered by the simulation model that must be noted. (See Appendix A for a depiction of 


the ECB used in the simulation.) These assumptions and limitations include: 


Transition probabilities to and from certain states within the state space 
diagram are based solely on the author’s professional judgment. Since the 
EFGOM 1s a proposed weapon system, detection, performance and 
survivability data are not currently available. 


Forward Observer (FO)/Sensor performance is not based on any technical 
data. No thermal or image intensification devices were modeled. It is only 
assumed that the FO/Sensor has a finite range (10 km) in which it can detect 
targets, and a defined sector of scan. 


The FO will attempt to establish communications with BDE or AFATDS to 
pass a CFF message only a limited number of times. After a defined number of 
attempts, he will drop the acquired target, and go back to searching and try to 
establish communications upon his next acquisition. This assumption ts based 
on professional experience that would indicate that a FO would not stop doing 
his mission (scanning the battlefield for targets) just because he could not 
contact his higher headquarters. 


High value targets are assumed to be traveling in march column to the location 
of U.S. forces at a uniform rate of march, regardless of any terrain 
considerations. Additionally, the vehicles are dispersed in the march column in 
exact doctrinal distances between vehicles, platoons, companies, battalions, 
regiments, etc. This march column ts based on the march column in HRS 33.5. 


Once a target has be acquired by a FO, the FO 1s able to track that target 
throughout the engagement, and will not CFF on that target again. 


All targets will enter the main battle area by traversing through a single 
EFOGM engagement area (EA). The five sensors are arranged around the EA 
such that three are placed on the leading edge of the EA and two are in depth 3 
km behind the leading edge of the EA. 


EFOGM batteries are assumed to be located so that the distance from the 
battery to the leading edge of the EA is the maximum range of the EFGOM. 
Additionally, all EFOGM fires are considered to have come from a single 
“super-battery”’, and not from tactically dispersed, individual gun systems 
controlled through a battery or platoon headquarters unit. 


The C’ node located at BDE or AFATDS is able to track the continuous 
location of a targeted high value target. This prevents multiple CFFs from 
being processed on the same target. Since the number of high value targets is 
limited within the array of all targets, the ability to track this location is 
feasible. 


While the digitized scenario depicted by CSC does not include a mission 
confirmation of the AFATD’s CFF, it is assumed that at a minimum, BDE will 
pass a mission confirmation and a “Fire/Hold Fire’ command to subordinate 
elements to be echoed down the chain of command. Without this 
consideration, a CFF feasibly could reach the battery, and not the BDE or PL, 
and be processed for firing by the battery. Since “silence equals consent” in 
the CSC model, the battery would assume that the lack of any 
acknowledgment from BDE meant to execute the mission, when in fact, it 
could be because BDE never received the mission. Considering the amount of 
attention focused on fratricide in today’s Army, it is illogical to think 
confirmation would not ever be required. 


Missile damage is not a function of target vehicle type, but purely a function of 
a Monte Carlo draw. 


Missile rates of fire are not considered. Since all guns are consolidated into 
one “super-battery” simultaneous or near simultaneous launches may occur. 


Shoot downs of EFOGM missiles while in flight, or EFOGMs that become 
erratic and lost in flight (LIF) are assumed to occur at a time value equal to 
one half of the total time of flight of the missile to the intended target. 


No re-supply occurs during the battle. The EFOGM battery fights with only 
its original basic load of missiles. Considering possible scenarios for early 
entry forces, this is a distinct possibility. 


For simulation purposes, the lead element 1s located 9 kilometers from the 
leading edge of the EFOGM EA at the start of the simulation. This allows for 
the earliest possible acquisition of targets in the simulation. 


Il. STATE SPACE DESCRIPTION 


A. THE EXTENDED CLOSE BATTLEFIELD 


To study the effects of the EFOGM on the Extended Close Battlefield (ECB), a 
description of events that occur in that area of the battlefield is necessary. A description 
of those battlefield events could be considered as a collection of random variables X(t), 
where X(t) is the state of the battlefield at time t. This description of the battlefield can be 
considered a stochastic process, where the state space of that stochastic process is defined 
as the set of all possible values that the random variable could assume. A definition of the 


stochastic process 1s provided in Jntroduction to Probability Models by Sheldon M. Ross. 


... consider a stochastic process {X,,n = 0, 1, 2,...} that takes on a finite 
or countable number of possible values. Unless otherwise mentioned, this 
set of possible values of the process will be denoted by the set of non- 
negative integers {0, 1, 2,...}. If X, =i, then the process is said to be in 
State i at time n. We suppose that whenever the process is in state 1, there 
is a fixed probability P,; that it will next be in state 7. That is we suppose 
that 


P{Xn+1= jl Xn=1, Xn=1 = eles | Sh 0 0p — (1) 


for all states ip, 4),..., in1, 7,7 and alln 20. Such a stochastic process is 
known as a Markov Chain. Equation (1.1) may be interpreted as stating 
that, for a Markov chain, the conditional distribution of any future state 


Xn+i given the past states Xo, X1, . . ., Xn: and the present state Xj, is 
independent of the past states and depends on the present state. [Ref. 10, 
pai35) 


Note that for a discrete time Markov chain, transition times are all one time step. This 
particular Markov chain refers to a process with stationary probabilities of transition (.e., 
the probability of transitioning from one state to another is not a function of the previous 
State, time, range or any other factor). Such a description of the ECB could be used, but 
it would not account for such factors as ranges to targets that are critical when 


considering missile engagements. Therefore, the classical Markovian approach can not be 


used since state space transition probabilities and sojourn times are functions of previous 
States and do not posses the memoryless property, therefore prohibiting a closed form 
solution. Likewise, a continuous time Markov Chain would be inappropriate since the 
conditional distribution of a future state given the current state and past states depends 
only on the present state, and 1s independent of the past. This would not be appropriate, 
for example, in a state that describes the processing of a call for fire (CFF). CFF 
processing 1s dependent on when the CFF was initiated, which is a past event. However, 
the notion of a Semi-Markov process can be used to generalize the ECB state space. 
Again from Ross: 

Suppose that a process can be in any one of N states 1, 2,..., N, and that 

each time it enters state 7 1t remains there for a random amount of time 


having mean ul; and then makes a transition into state 7 with probability Pj. 
Such a process 1s called a Semi-Markov process. [Ref. 20, p. 325] 


For the state space description of the ECB, the Semi-Markov notions of the 
probability of transitioning from one state to another (Pj;'s), the transition (sojourn) times 
from State to state (Tj;’s), absorbing states, total number of visits to specific states, and the 
proportion of time that the model spends in a particular state i (P;) will be used to develop 


measures of effectiveness within the state space description. 


B. STATE SPACE INTERPRETATION 


1. General 


The state space used in this thesis replicates the area in the ECB where FOs are 
positioned to provide targeting information to the controlling DRB headquarters. For the 
purposes of this thesis, an FO represents any and all types of information gathering 
systems from ground emplaced sensors, radars, airborne platforms and scout vehicles to 
dismounted forward observers. It consists of j + 24 states, where j is defined as the 
number of high value targets determined to be present in an attacking enemy force. What 


follows 1s a description of each state, the explanation of transition probabilities from that 
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state, and how sojourn times are determined for that state. A graphical depiction of the 


state space 1s shown 1n Figure 1. 
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Figure 1. State Space Diagram 
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2 Sensor States 


States One through j describe the FO scanning the battlefield, acquiring targets, 
and the possible detection and destruction of the FO by the enemy. The description of 


each state Is: 


e State One - FO Initialization. In State One the FO 1s “initialized” or brought to 
life. Each FO has a maximum visibility range and scans a sector based upon 
user defined scan arcs and fields of view (FOV) based upon the notion of the 
Cartioid Distribution [Ref. 4, p. 80]. The observer/FO may transition to either 
State Two (FO Detected) or State Three (FO Undetected / Not Destroyed). It 
is assumed there is no time needed for transition from initialization to 
searching, and that the chance of detection by enemy forces remains constant 
throughout the simulation. 


e State Two - FO Detected. In State Two, enemy forces have detected the FO. 
The FO may transition from State Two to State Three (FO Destroyed - an 
absorbing state) or State Four (FO DetectedlEngaged but Not Destroyed). It 
is assumed that once detection is made, the FO completes the CFF and is then 
destroyed. 


e State Three - FO Destroyed. State Three is an absorbing state that describes 
the detection of the FO by enemy forces and its subsequent destruction. There 
is no transition from State Three. 


e State Four - FO DetectedlEngaged but Not Destroyed. State Four is a transit 
node from State Two (FO Detected) to State Six (FO Searches). The 


transition probability is 1.0 (if not destroyed the FO will search), and the 
transition time is 0.0. 


e State Five - FO Undetected and Not Destroyed. State Five is a transit node 
from State One (FO Initialization) to State Six (FO Searches). The transition 


probability is 1.0 (if undetected and not destroyed the FO will search), and the 
transition time is 0.0. 


e State Six - FO Searches. In State Six the FO begins its search of the battlefield 
for high value targets within its FOV. Should it detect a high value target, it 
will transition to the specific detection state (7 through /) with probability 1.0. 
The sojourn time from the search state is a function of a detection algorithm. 
However, if the FO fails to detect a target within a time limit, it transitions 
back to State One with a transition time equal to the time limit. This transition 
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possibility provides for a continuous probability that the FO could be detected 
throughout the simulated battle. This is a realistic depiction of the searching 
that occurs on the battlefield by both friendly and enemy forces. 


States 7 through 7 - Detect High Value Target. States 7 through j describe the 
detection of a specific high value target by an FO. Each state will transition to 


State 7 + | (Transmit a Call for Fire) with zero sojourn time and probability 1.0 
if detection occurs for that specific target. Once a target is detected, it is 
assumed the FO will immediately call for fire (CFF) on the detected target 


States j + 1 through j + 11 describe the action of the FO transmitting a CFF to the 
DRB TOC to launch of the EFOGM: 


Sy. 


State 7 + 1 - Transmit Call for Fire (CFF). State 7 + 1 describes the FO 


transmitting a CFF on the detected target. There 1s a chance that his CFF will 
not be received at the DRB HQ or the AFATDS (State 7 + 2) or that the DRB 
HQ/AFATDS will receive the CFF and begin processing the mission (State j + 
3). The sojourn time in state 7 + 1 1s negligible. Either radio contact is made 
immediately, and transition is made to state j + 3, or the FO fails to make 
contact and transition 1s made to state j + 2. 


State j + 2 - CFF Not Received/Re-Transmit CFF. Should the FO fail to make 
contact with the DRB HQ/AFATDS, he will continue to attempt 
communications for a period of time before abandoning the attempt, returning 
to search for other targets (a return to State One). It is assumed that each 
attempt to re-contact the DRB HQ/AFATDS will take approximately five 
seconds. Should contact be made, the process continues to State j + 3 (CFF 
Received/Processing Begins). 


State 7+ 3 - CFF Received and Processed at DRB S2/FSO Cell or at 
AFATDS. State j + 3 accounts for the largest portion of time consumed in the 


state diagram. The probability of transitioning from this state to State 7 + 4 
(EFOGM Launch) depends on the time required to process the CFF from the 
he: 


Processing States 


A specific discussion of the time required to process the CFF depends on the 


scenario that is being replicated. Specifically, there is a different logical process for the 


two types of C° systems. Processing for both the digitized and non-digitized systems will 


now be discussed. (Since the user decides to execute either a non-digitized or digitized 


scenario, there is a duplication of states / + 3a through j + 37.) 
a. Non-Digitized Processing 


States / + 3a through / + 3f depict the non-digitized processing of the 
CFF. These states depict the expanded version of the single state +3 shown in Figure 1. 


The states are depicted in Figure 2, and will now be discussed in detail. 
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Figure 2. Non-Digitized Processing Sub-State Space 





e State 7 + Ja - CFF Received and Processed at DRB $2/FSO Cell. Time begins 
when the FO transmits the first CFF to the DRB HQ. There a radio telephone 
operator (RTO) must stop or finish any job he is currently working on, record 
the CFF, and decide what action is needed. The target must then be passed to 
the brigade fire support element (FSE). This process takes a varying amount 
of tme, depending on the nature of the target. If either the brigade fire 
support officer (FSO), or fire support non-commissioned officer (FSNCO) is 
not busy with other fire missions, then they will process the mission 
immediately. If busy, the target data are “queued”, and will be rechecked when 
either is free and all previously “queued” target information is processed. The 
FSO, or FSNCO then decides if the target is a valid EFOGM target according 
to the brigade commander’s intent. Once the target 1s processed, it is delivered 
to the EFOGM liaison officer (LNO) for processing. If the LNO is not 
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immediately available, again the target data are queued and will be checked 
when the LNO is free. 

Should the total processing time from the beginning of the FO’s CFF to the 
beginning of the LNO’s processing exceed a maximum allowed time, the target 
may be dropped. For example, targets moving at 5.55 meters per second can 
travel over two kilometers in six minutes. Since additional processing time will 
continue to accumulate after the target 1s sent to the LNO, this may seriously 
decrease the EFOGM probability of target acquisition in the engagement area. 
Should this occur, States 7 + 3a or 7 + 3b could transition into State j + 3c 
(Excess Process Time), which ts an absorbing state. If not, the EFOGM LNO 
verifies the mission, determines which EFOGM asset will fire the mission, 
prepares the launch order and attempts to pass the fire mission by FM radio to 
the EFOGM battery (State 7 + 3e). In reality, the LNO would also determine 
the number of launchers and number of missiles to be fired for the mission. 
However, for purposes of this simulation, it is assumed that only one missile 
will be fired per target, and there 1s only one consolidated battery containing all 
missiles. 


State j + 3b - Transmit Fire Mission to EFOGM Battery. Similar to State j + 
1, this state represents the LNO attempting to pass the fire mission to the 


EFOGM battery over FM voice radio. The transmission may be received 
(State j + 3e), or not received and must be tried again (State 7 + 3d). Time 
Accumulation ts the same as State j + I. 


State j + 3c - Excess Process Time. State j + 3c is an absorbing state 
describing a CFF that has required too much time to process, and 1s dropped 
from consideration. 


State j + 3d - CFF Not Received/Re-transmit CFF. Same as State j + 2, 
however, should total processing time exceed a set standard, the target will be 


dropped and a transition into State j + 5 occurs. 


State 7 + 3e - CFF Received at EFOGM Battery/Battery Processing. In State j 
+ 3e, the battery platoon leader has received the fire mission from the LNO, 


and begins processing the fire mission. The platoon leader determines the aim 
point and missile way point, based on the terrain, target range, and target rate 
of march. Once the platoon leader projects an intercept zone, he prepares the 
fire mission and issues the launch order to the firing squads. 


State j + 3f - Battery Applies Data to Launcher. As data arrive at the launcher, 


the gunner inputs the data into the launch computer and executes the launch 


— 
tn 


command. The cumulative processing time becomes the sojourn time to State 
J +4 (EFOGM Launch). 


b. Digitized Processing 


States j + 3a through j + 3/ depict the digitized processing of a CFF. 
These states depict the expanded version of the single state +3 shown in Figure 1. The 


states are shown in Figure 3, and wil! now be discussed in detail. 


e State 7 + 3a- CFF Received and Processed by AFATDS. Once the target is 
received by AFATDS, it compares the target to the target priority, assets 
available to fire and the brigade commander’s plan. It computes which units 
should fire, and the number of missiles for the mission. It then simultaneously 
passes the CFF to the BDE FSE, the platoon leader and the battery. 
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Figure 3. Digitized Processing Sub-State Space 


e States 7 + 3b/d/f - CFF not received by BDE/PL/Battery - AFATDS re- 
transmits. In these three states the CFF transmitted by AFATDS was not 
received by some or all of the receiving units. AFATDS will continue to 
attempt the transmission until contact 1s made. 
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State j + 3c - AFATDS mission received by brigade. Once the mission arrives 
at the brigade level, it is checked to ensure that an acceptable firing solution 


was reached by AFATDS. While the CSC report states that “silence equals 
consent’, it is felt that such a hands-off approach to engagements is unlikely, 
especially with the growing fratricide concerns within today’s armed forces. 
Therefore, in this state space depiction, brigade confirms the mission and 
passes a “Fire” or “Hold-Fire” to the platoon leader. 


State 7 + 3e - AFATDS mission received by PL. The platoon leader has 
received the mission specifics from the AFATDS and awaits the command 


from BDE to hold fire or fire the mission. 


State 7 + 3g - AFATDS mission received by battery. Once the mission arrives 
from AFATDS, the battery immediately begins processing the data into a firing 


solution. The battery waits until the fire command is echoed by the platoon 
leader to the battery to apply the firing data to the launcher. 


State 7 + 3h - Mission confirmation not received by PL. Similar to States j + 
3b/d/f, BDE will continue attempting the transmission until contact is made 


with the PL. 


State j + 31 - PL receives mission confirmation from BDE. Once the PL 
receives the mission confirmation from BDE, he immediately passes the “Fire” 


or “Hold-Fire’” command to the battery. 


State 7 + 3) - Mission confirmation not received by Battery. Similar to State j 
+ 3h, the PL will continue attempting the transmission until contact is made 


with the battery. 


State j + 3k - Battery completes mission data processing. In this state the 
battery has finished transferring mission data to the system and selecting and 


validating the missile’s route to the target. The battery waits here for mission 
confirmation from the PL. 


State j + 3/ - Battery applies data to the launcher. Once the PL’s confirmation 
arrives, the gun applies the data to the launcher, and executes the launch 


command. 
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3. Missile Launch States 


States ] + 4 through j + & describe the possible occurrences at missile launch, and 


are described as follows:. 


e State 7 + 4-EFOGM Launch. State j + 4 describes the gunner executing the 
launch command (1.e., “pulling the trigger’). State 7 + 4 can transition into one 
of three states: 7 + 8 (Successful Launch), 7 + 6 (Misfire/Re-shoot) or j + 5 
(EFOGM Lost). There is no time involved with any of these transitions. 


e State 7 + 5 - EFOGM Lost. This is an absorbing state describing a missile that 
leaves the launcher and fails to follow launch commands (e.g., an erratic 
round). 


e State j + 6 - Misfire/Re-shoot. State j + 6 describes a missile that fails to leave 
the launcher once the launch command has been executed. There are misfire 
procedures which the gunner must follow, and these account for the transition 
times from this state to four other possible states: 7 + 5 (EFOGM Lost), 7 +7 
(Dud Missile/New Missile selected), 7 + 8 (Successful Launch), or the state can 
transition back into itself. 


e State 7 + 7 - Dud Missile/New Missile Selected. Should a missile continue to 
misfire, it may in fact be a dud round that will not fire. Should this occur, a 
new missile will be selected for launch and the gunner will “‘pull the trigger” for 
that missile. The time needed to select a new missile and apply launch data is 
the transition time from this state to states 7+4,7+5orj+8. Again, 
transition probabilities are defined by the user. 


e State j + 8 - Successful Launch. This state describes the missile that has left 
the launcher and is flying properly. This state concludes the launch states and 
transitions to the flight states of 7 + 9 (EFOGM Detected in Flight/Engaged). 
and j7 + 10 (EFOGM Undetected in Flight). The transition time for missile 
flight time is not considered from this state. 


4. Missile Flight and Impact 


States j + 9 through j + 16, describe possible events during the flight of missile, 
impact of the missile on the target, or failure of the missile to impact on the target, and 


have the following criteria: 
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State 7 + 9 - EFOGM Detected and Engaged. Within this state the EFOGM 
missile has been detected and engaged by enemy forces. Destruction of the 
missile causes a transition into State 7 + 11 (EFOGM Destroyed), while failure 
of the enemy to destroy the missile transitions the process into State 7 + 12 
(EFOGM Missed by Enemy). The only calculated transition time would be 
from j + 9 to 7 + 11, that is estimated to occur at one half of the missile’s flight 
time to the target. 


State j + 10 - EFOGM Undetected in Flight. This state indicates that the 
missile has not been detected in route to its target. However, it could become 
erratic during flight and be “lost” (State 7 + 13). If this does not occur, the 
state space transitions to State j + 14 (EFOGM Target Acquisition/Lock On). 
No transition times between states are considered here unless the missile is lost 
in flight, which is estimated to occur at one half of the flight time to the target. 


State j + 11 - EFOGM Destroyed. This is an absorbing state describing the 
destruction of the EFOGM by enemy forces. 


State 7 + 12 - EFOGM Missed by Enemy. This state represents an enemy 
failure to destroy the EFOGM in flight. The missile may now acquire and lock 


onto an enemy target (State 7 + 13), or become lost in flight (State 7 + 12). 
Transition time is calculated only if the missile becomes lost in flight. This 
transition time is estimated as one half the time of flight to the target. 


State 7 + 13 - EFOGM Lost in Flight. This is an absorbing state describing a 
missile, in flight to a target, that becomes erratic or has a malfunction that 
causes it to fly improperly, and become lost. This is assumed to occur at a 
time equal to one half the time of flight of the missile to the intended target. 


State 7 + 14 - EFOGM Target Acquisition/Lock On. At this point, the 
EFOGM has “tipped over” so that its seeker head is looking into the target 


area and has acquired/locked on to a target vehicle. The missile can either 
impact on the target (State j + 15), or miss the target (State 7 + 18). Transition 
time from State j + 14 to either 7 + 15, or 7 + 18 1s equal to the time of flight of 
the missile to the target. 


State 7 + 15 - EFOGM Impact. In this state, the missile has impacted on the 
target vehicle. There are two possible transitions from this state: either the 
missile causes catastrophic destruction of the target (State 7 + 18), or the target 
is not destroyed upon impact (State 7+ 17). There is no transition time from 
this state to either state 7 + 17 or j + 18. 
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State 7 + 16 - EFOGM Misses the Target. The missile has failed to impact the 
target in this state. However, it is assumed that despite a failure to impact on 
the target, there is a possibility of some collateral damage due to the explosion 
of the warhead if the missile lands close enough to the target. Transition time 
from this state to one of the States 7 + 19 through j + 22, is considered to be 
zero. States 7 + 19 through j + 22 describe the levels of damage that could be 
caused by a ‘near miss’ of the EFOGM. It is also assumed that a ‘near miss’ 
can not cause catastrophic target destruction (State j + 18). 


State 7+ 17 - Target Not Destroyed. State 7 + 17 indicates that though the 


EFOGM impacted on the target, the target was not destroyed. State 7 + 17 
may transition to any one of the States j + 19 through j + 22 that describe the 
varying levels of damage caused by the missile’s impact. There is no transition 
time between these states since the time between impact and a specific level of 
damage occurring is considered to be zero. 


Battle Damage and Damage Assessment 


States 7 + 18 through j + 22 identify the level of damage, including destruction, 


caused by the EFOGM on the target, while States j + 23 and j + 24 describe the battle 


damage assessment (BDA) made by the FO on the target: 


State j + 18 - Target Destroyed. State j + 18 is accessible only from State j + 
15 (EFOGM Target Hit). It can transition to one of the BDA states, 7 + 23 
(Re-shoot Mission), or j + 24 (Go to New Target). It is assumed that even 
though the target is catastrophically destroyed, the FO may make an incorrect 
BDA of the target and execute a repeat CFF on that target (a transition to 
State j + 23). 


States 7 + 19 through 7 + 22 - Damage Level States. These states describe the 
level of damage caused by the impact, or near miss of the EFOGM. There are 


five possible damage levels: Destroyed, Heavy, Moderate, Light, and None. 
The transition probabilities from either State j + 16, or j + 17 to these states 
will vary, since it is assumed, for example, that the probability of Heavy 
Damage (State j + 19) would be much higher from State 7 + 17 (Target Hit but 
Not Destroyed) than from State j + 16 (Target Missed). All damage level 
states transition to States 7 + 23 (Re-shoot Target) or j + 24 (Go to New 
Target). 


State 7 + 23 (Re-shoot Target). In this state, the FO has determined that the 
level of damage is insufficient enough to warrant a repeat of the fire mission. 
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The transition time from the damage states (/ + 16 andj + 19 through j + 22) 
to this state represents the time needed by the FO to asses target damage and 
choose a course of action. Should a transition to this state occur, the only 
transition out 1s to State 7 + 1 (Transmit CFF) with probability of 1.0. 


e State 7+ 24 (Goto New Target). In this state, the FO has determined the 


damage level 1s sufficient not to re-fire the target. The FO then begins 
searching for new targets (State 1). The transition time to this state 1s defined 
as the time needed by the FO to asses target damage and choose a course of 
action. 


c. USE OF THE STATE SPACE DIAGRAM 


As mentioned, the state space diagram 1s used to model a portion of the Extended 
Close Battlefield. The issue then becomes how to simulate this state space to return 
values of interest; specifically, the times required to process CFFs and prosecute enemy 
targets. The general layout of the state space lends itself for use in a simulation model that 
can replicate events on the battlefield and the times involved with those events. The 
following chapter will discuss the nature of that model, and the methods used for 


extracting desired information from the state space depiction. 
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Il. MODEL DESCRIPTION 


A. GENERAL DESCRIPTION 


The simulation model used for this thesis was written in PASCAL. It uses a 
uniform probability random number generator to replicate a Monte Carlo draw which is 
used to determine along which arc the state space simulation will advance. A description 
of the random number generator is contained in the next section. Monte Carlo 
comparisons are based upon user defined transition probabilities (P;’s). Additionally, 
normal and lognormal probability distribution random number generators are used to 
determine the transition times of several states within the simulation. For the purposes of 
this thesis, certain transition probabilities and times are based on the author’s professional 
judgment. A complete listing of transition probabilities and their source ts located in 
Appendix E. However, the majority of the transition times are based upon previous 
studies that define specific actions and functions of the command and control process as 
time distributions [Ref. 2, p. 2-14]. PASCAL was chosen as the simulation language due 
to its ease of use, and because it does not require a particular operating system, other than 
DOS, to function. This ensures its execution on stand alone personal computers present in 
the majority of U.S. Army analytic agencies today. 

Output from this simulation will allow the user to collect the number of times 
certain states were “visited” and the times from target detection to target prosecution and 
missile interaction. With data from a large number of simulation runs of both the digitized 
and non-digitized C° systems, operating range values for the two systems can be obtained. 
Summaries of these data can then provide possibie input parameters for the EFOGM 
system in other simulation models, or can be used independently to study the differences in 
performance of the digitized versus the non-digitized system altering various input 


parameters. 
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B. RANDOM NUMBER GENERATOR 


The computer code shown in Appendix C is the method used within the simulation 
model to generate random numbers from the uniform distribution on the interval [0,1] 
(noted as U(0,1)). The importance of the U(0,1) generator is that the random variates 
from other probability distributions (normal, lognormal, beta) can be generated by 
transforming a U(O,1) variate in a manner determined by the desired distribution [Ref. 7, 
p. 421]. The method used to generate the normal and lognormal variates is the Inverse 
Transform technique. An example from Law and Kelton provides a simple explanation of 


the technique: 


Let X have the exponential distribution with mean B . . . The distribution 
function 1s: 


rel ee 310) 
Foe). ee 
0 otherwise 


so to find F’, we set u = F(x) and solve for x to obtain: 
(7) = ole 


Thus, to generate the desired random variate we first generate a U~U(0,1) 
and then let X = -B In U. [It is possible in this case to use U instead of 1 - U, 
since 1 - U and U have the same U(0,1) distribution .. .] [Ref. 7, pg. 466] 


Transition time values from the normal and lognormal distributions were 
generated, in this manner, for use in the simulation. To ensure that the random number 
generators were functioning properly, a small test program was written to generate 1000 
random variates from the normal distribution on the interval (0.25, 0.2778] and from the 
lognormal distribution on the interval [0.41667, 0.0611]. These two values are provided 
in the BEWSS model as the random times to transmit a voice radio message and FSO 
processing time (in minutes), respectively. Shown in Appendix D are graphical and 
analytical distribution fits for each sample of 1000 variates. These were produced by the 


Statistical program A Graphical Statistical System (AGSS). As shown in the Appendix D, 
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the samples follow their respective distributions very closely, and verify the validity of the 


random number generators. 


c SIMULATED FLOW OF THE BATTLE 


1. Main Program 


The main program is a simple loop that executes a user defined number of 
replications of either the digitized, or non-digitized simulation model. Within this loop, a 
“master clock” is initialized that tracks current “battle time” throughout the program. A 
queue of records stored as a linked list (an event queue), called the Calendar, is created 
which stores Events to be executed, sorted by execution Time. Parameters such as which 
observer is searching (Obsvr), the target number acquired and engaged (TgtNumb), the 
Missile engaging a target, the Launcher from which the missile was fired, and a pointer 
to the next record (Next), are also stored within each record of the Calendar. For non- 
digitized replications, additional queues for CFFs that the FSO and LNO must process 
when they become free from processing previous CFFs, are initialized. The first event, 
‘“StartBattle” is placed on the Calendar to be executed at time 0.0 using SCHEDULE, 
which 1s a procedure outside the main program, but warrants explanation now. 

The SCHEDULE procedure controls inputs to all queues. The first check 
SCHEDULE makes is to ensure the event time about to be scheduled is greater than, or 
equal to the current clock time. (Obviously, an event can not be scheduled to occur in the 
future if it has an execution time in the past.) If the event time meets the logic check, the 
event 1s placed on the queue. For example, scheduling the Event “CFF” (Call For Fire) 
requires values for the parameters Observer and Target (which are both integer values 
from one to the maximum number of FOs, or maximum number of targets, respectively). 
Because the Missile and Launcher parameters are not needed in the event “CFF”, they 
receive integer values of zero (0). The parameter EventTime is a real value representing 


battle time in minutes. For the Calendar, FSOList, and LNOList, SCHEDULE places 


the values for the parameters listed above on the queue in a time sorted manner. The 
model includes calls to SCHEDULE when it determines an event should occur. 

Finally, the main program calls either MASTERTIMER or MASTERTIMER2. 
MASTERTIMER is the subprogram for a non-digitized replication, while 
MASTERTIMER2 replicates the digitized C’ system. These two sub-programs remove 
and schedule items from and to Calendar, FSOList and LNOList, keep track of the 
simulation clock, and execute the simulation. Each has several subroutines. All 
subroutines, minus those that replicate the digitized, or non-digitized form of processing, 
are similar for both sub-programs. Discussion of these subroutines will cover the common 
initial subroutines, break out the specifics for the digitized and non-digitized cases, and 
then discuss the final common subroutines. For simplicity purposes, future references to 
subroutines or procedures that are peculiar to either MASTERTIMER or 
MASTERTIMER? will be noted as (ND: non-digitized) and (D: digitized), respectively. 
Additionally, some procedures will use probabilistic time distributions from which various 
processing and wait times for the simulation are drawn. These times are taken from the 


CSC technical report [Ref. 2], and are summarized in Appendix E. 


28 Initial Common Subroutines 


As stated, MASTERTIMER and MASTERTIMER2?2 control program execution. 
While the Calendar is not empty, they remove the parameters from the first record in the 
queue using the procedure “DeQueue’”’. They then compare the event name taken off the 
Calendar to a list of indexed sub-routine event names. Once a match is found, the sub- 
program takes the index number for that event name and searches its indexed subroutines 
until a match is found between the index number and the indexed sub-routine. It then 
executes the events in the matching indexed sub-routine at the time drawn from the 


Calendar. The specifics of the initial common subroutines will now be discussed. 


26 


a. “StartBattle”’ 


This sub-routine initializes any globai variables to their beginning of battle 
values. For the computer model used in this thesis, the user defines these starting 
parameters by altering input parameters located in a text data file (example shown in 
Appendix F). “StartBattle’’ reads in the user defined parameters from the data text file, 
and as it loops over every combat entity, it is initialized to its beginning battle value. It 
also assigns transition probabilities, and inputs the initial ranges from every high value 
target in the scenario to every FO and EFOGM launcher. In particular, it schedules the 


sub-routine ‘Initialize’ for all FOs to start at battle time 0.0. 
b. “Initialize” 


This sub-routine corresponds to States One through Five in the state space 
diagram (see Figure 1). It samples a U(O,1) and compares the sample to indexed 
transition probabilities that were assigned in the “StartBattle” sub-routine, to determine if 
the FO is detected. If not detected, it schedules the sub-routine “Search”. If detected, the 
sub-routine makes another U(0,1) draw to determine if the FO 1s destroyed. If not 
destroyed, it schedules a “Search” and returns control to MASTERTIMER or 
MASTERTIMER2. If destroyed, it decrements the number of remaining FOs, updates the 
Status of that particular FO to “Destroyed” and returns control to MASTERTIMER or 
MASTERTIMER?2. 


CG “Search’’ 


“Search” corresponds to States 6 through j in the state space model. For 
the searching FO, it loops over all possible targets. If the damage level to the target is not 
“Destroyed”’, it draws a U(O,i) and checks the current range between the FO and that 
target. If the U(O,1) draw results in the target being in the FO’s field of view (FOV) and 
the range to the target is within the maximum range of the FO, but not behind the FO (at a 


negative range), and the FO has not previously transmitted a CFF for this particular target, 


27 


a time to acquire the target is determined using the DYNTACS detection rate equation 
[Ref. 4, pp. 77-79]. (See Appendix G for a discussion of the DYNTACS equation and 
related detection algorithms.) If that acquisition time is the smallest over all acquired 
targets for that FO, it is stored, along with the corresponding target number. If not, it 
loops to another target. This continues until all targets are sampled. Once the loop over 
all targets 1s complete, the minimum acquisition time for that FO over all targets is 
compared to the maximum time any FO can search without being re-initialized. This 
check prevents an FO from searching for long periods of time without running the risk of 
re-detection. If the minimum acquisition time is less than the maximum cycle time, a 
“CFF” (Call for Fire) is scheduled at the clock time plus the calculated minimum 
acquisition time. If not, the FO is re-initialized at the clock time plus the maximum cycle 


length time, and control returns to MASTERTIMER or MASTERTIMER2. 


d. “CERF” 


This sub-routine covers states j +1 andj +2. A U(0,1) is drawn to 
determine the success of the attempted contact between the FO and the DRB HQs (ND) 
or AFATDS (D). If successful, a “BdeProcess” (ND) or “AFATDS” (D) is scheduled at 
the input clock value plus a random time value from the normal distribution, and control 
returns to MASTERTIMER or MASTERTIMER2. If a U(0,1) draw indicates a failure to 
successfully contact the DRB HQ/AFATDS, additional attempts are made using different 
U(O,1) draws. This process continues until the Monte Carlo draw indicates successful 
contact is made, or until a user defined amount of time has past. If that point is reached, 
the FO quits trying to make contact and an “Initialize” is scheduled at the input clock 
value plus the user defined time, and control then returns to MASTERTIMER or 
MASTERTIMER?2. 
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3; Non-Digitized Processing Common Subroutines 


a. “BdeProcess”’ 


“BdeProcess” is one component of State 7 + 3 (see Figure 2). Here, a 
random lognormal value is drawn to represent the ume necessary for an RTO to process 
(record) the FO’s CFF. A truncated random normal value is then determined to represent 
the time spent checking the status of the FSO. These two events constitute the brigade 
processing time. The event “FSOProcess” is then scheduled on the Calendar at the clock 
value plus the brigade processing time. If the FSO can process no more CFFs (1.e., there 
are no missiles left to shoot), the message 1s not passed to the FSO. Control] then returns 


to MASTERTIMER. 
b. “FSOProcess” and “FSOFree”’ 


Another component within State 7 + 3 1s the “FSOProcess” sub-routine. If 
the FSO is busy when the event is taken off the queue (checked with a boolean), the 
mission is queued on the FSOList. If the FSO is not busy when the CFF arrives, two 
random lognormal samples drawn to replicate the times needed for the FSO to process the 
CFF and then deliver it to the EFOGM LNO cell. If the total process time for the fire 
mission, from the FO’s CFF to hand-off to the LNO, is less than a user defined value, the 
mission is delivered to the LNO by scheduling “LNOProcess”’. Should the total process 
time exceed the user defined amount, the CFF is not processed any further. It is assumed 
that because of the delay in processing that the additional time required to process the 
mission and then fly the missile to the target will seriously degrade the possibility that the 
EFOGM will acquire the proper target within the engagement area. At the same time the 
“LNOProcess” is scheduled, an “FSOFree”’ is scheduled to occur. The ‘““FSOFree” 
procedure changes the boolean value for the FSO and checks the FSO queue for any 


missions. If any missions are on the queue, they are taken off and an “FSOProcess”’ is 


ie 


scheduled at the current clock value. If the FSOList is empty, the FSO waits for his next 
CFF. Control then returns to MASTERTIMER. 


a “LNOProcess”’ and “LNOFree”’ 


The “LNOProcess” sub-routine first checks to ensure that the LNO is not 
busy processing a previous CFF. If he is, the new mission is placed on the LNOList and 
control returns to MASTERTIMER. If the LNO is not busy and the total number of CFF 
processed is less than the maximum number of missiles, two random lognormal samples 
are drawn representing the EFOGM LNO processing time and the time needed to transmit 
the firing order over the radio to the battery. A “LNOFree” is then scheduled at the sum 
of the two samples. 

Similar to the “FSOFree’’, the “LNOFree” procedure changes the boolean 
value for the LNO and checks the LNO queue for any missions. If any missions are on the 
queue, they are taken off and an “LNOProcess” is scheduled at the current clock value. If 
the LNOList is empty, the LNO waits for his next CFF 

A Monte Carlo draw is then conducted to confirm if the EFOGM battery 
can receive the CFF from the EFOGM LNO by radio. If contact can be made, a 
“BtryProcess” 1s then scheduled at the clock time plus the sum of the two lognormal 
samples, and control returns to MASTERTIMER. Should the Monte Carlo draw indicate 
that radio contact is not initially made, U(O,1) draws are made until the value of the draw 
is equal to, or exceeds the value needed to transition into State 7 + 3e (CFF Received by 
EFOGM Battery). Each retransmission adds a user defined amount of time to the total 
processing time. Once contact is established, a “BtryProcess” is scheduled at the sum of 
the two lognormal samples and the extra time needed to establish radio contact. Control 


then returns to MASTERTIMER. 
d. “BtryProcess” 


The “BtryProcess” sub-routine compiles the total time for the battery 


platoon leader to predict an intercept location for the target, process and issue the fire 
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command, and for the EFOGM gunner to apply the data to the launcher. These four times 
are represented by four distinct lognormal random number draws. By calling the 
procedure “GetLaunchTime’”’, the launch time for the mission is determined. If the target 
is currently within the maximum range of the missile, a “Launch” is scheduled 
immediately. If the range to the target, after the time incurred during a flight time to 
maximum missile range, is greater than the maximum range of the missile, the procedure 
calculates the time at which the missile should be launched, taking into account the flight 
time of the missile to maximum range, and the speed of the target vehicle to the 
engagement area. The event “Launch” is then scheduled for a time equal to the current 


clock value plus the time to launch. Control then returns to MASTERTIMER. 


4. Common Digitized Processing Subroutines 


a. “AFATDS” 


‘““AFATDS” represents the processing of a CFF by the AFATDS computer. 
First, the procedure “CFFTimes” calculates the amount of time (in terms of the number of 
transmission attempts needed to establish radio contact) it takes the AFATDS to establish 
radio communication with the BDE, PL and battery. It does this by sampling a U(O,1) and 
comparing it to the user defined probability of establishing first time contact. If the Monte 
Carlo draw is successful, “CFFTimes” returns the integer value “1”. If unsuccessful, 
U(O,1) draws are repeated until the draw is greater than the user defined probability of not 
establishing contact on succeeding tries. Each successive try increments a counter until 
contact is established. “CFFTimes” then returns the integer value of the number of 
attempts needed to establish radio contact. ““AFATDS” then checks to ensure that a CFF 
has not been previously processed for the particular target and that the total number of 
CFFs processed is less than the maximum number of EFOGM. If these conditions are not 
met, the sub-routine does not schedule the target for further processing. If the conditions 
are met, the total number of CFF is incremented, and the target is registered as being 


processed. The value of “CFFTimes” is sampled each time for the brigade, platoon leader 
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and battery. If the °CFFTimes” value is “1”, a “BdeProcess”’, “PLProcess” or 
‘““GunProcess’”’ is scheduled at the clock value plus the time required for the AFATDS to 
process the CFF, and then send it to a subordinate unit. Should the “CFFTimes” not equal 
‘*]”, it schedules a “BdeProcess’’, “PLProcess” or “GunProcess”’ at the clock value plus 
the sum of AFATDS processing time, transmission time, and the “CFFTimes” value 


multiplied by a user defined amount of time required for retransmission attempts. 


b. “BdeProcess” 


“BDEProcess” first samples the number of times needed to establish radio 
contact with the PL using the procedure “CFFTimes’’, discussed above. It then checks to 
ensure the acquired target is “High Priority”. If so, and the “CFFTimes” value is “1”, it 
schedules the subroutines “BDE1” and “PLProcess”’ at the clock value plus the time 
required for BDE to confirm the mission and transmit the confirmation to the PL. If the 
‘““CFFTimes” is not greater than “1”, the two subroutines are scheduled at the time just 
mentioned, plus the extra time required for retransmissions, as discussed in the section 
above. Should the target not be “High Priority’, “BDEProcess” schedules a “BDE2” and 


‘“PLProcess” using the same logic mentioned for “BDE1” and PLProcess”. 


Cc “PLProcess”’ 


Logically complicated, “PLProcess” accounts for the PL receiving data not 
only from AFATDS, but also mission confirmations from BDE. The PL normally receives 
initial mission data from AFATDS. However, if communications between AFATDS and 
the PL break down, it is assumed that the PL will accept mission information from the first 
unit, AFATDS or BDE, who sends him the new data. Therefore, there are three 
possibilities that could occur within “PLProcess’’. The first, and most likely possibility is 
that AFATDS passes the mission data to both BDE and PL simultaneously. “AFATDS” 
replicates this by scheduling “BDEProcess” and “PLProcess”’ to occur at the same clock 
value. Once the BDE and PL have received the mission data, ““BDEProcess” schedules a 


second “PLProcess” to occur at a later time, while the PL waits for BDE to confirm the 


mission. At the beginning of the second pass through “PLProcess”’ (scheduled by 
‘“BDEProcess’’), BDE has now confirmed the mission, and passes this confirmation to the 
PL. “PLProcess” simulates the PL confirming the mission, and then schedules a “PL1” Gf 
BDE passes a “Hold-Fire’’) or “PL2” (if BDE passes a “Fire’”’), and a “GunProcess’’. 
These subroutines are scheduled to occur at the clock value plus the time needed by the 
PL to contact the battery with mission confirmation and firing status. (Again, taking into 
account the “CFFTimes” number discussed previously.) 

The second possibility is that the PL does not receive the CFF data from 


AFATDS until after BDE passes mission information to the PL. “AFATDS” simulates 





this by scheduling a “PLProcess”’ at a time later than it schedules a “BDEProcess”. As 
“BDEProcess” is executed, it schedules a “PLProcess”’ at a time that is still less than the 
first “PLProcess scheduled by “AFATDS” . Should this occur, BDE has now initiated 
the CFF to the PL, and the PL will pass the mission confirmation and firing status to the 
battery, as mentioned in the previous paragraph. Finally, when the AFATDS mission 
arrives to the PL (simulated when the “PLProcess” scheduled by “AFATDS” is taken off 
the queue), the data are disregarded since it has already been forwarded to the battery. 
The final possibility occurs if the BDE does not receive the AFATDS data 
until well after the PL does. This case is no different, in a computer logic sense, than the 
first. The PL receives the data from AFATDS, and then must wait for BDE to confirm 


the mission, so he can pass the mission confirmation to the battery. 


a “GunProcess”’ 


‘“GunProcess” is also logically complicated, and consists of two possible 
cases, with two sub-cases within the first case. The first possibility is that the AFATDS 
mission data arrive at the PL and battery at the same time. “AFATDS” simulates this by 
scheduling the subroutines “PLProcess” and ““GunProcess” at same clock value. Battery 


processing of the data is replicated in “GunProcess” by sampling a series of three 
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lognormal time distributions, and changing the target number indexed boolean variable, 
ProcessingData, to “True”. There are two possible sub-cases for the first possibility. 

In sub-case one, the PL confirmation before the battery completes initial 
processing of the mission data, “GunProcess” samples the first two lognormal time 
distributions, and schedules a ““Gun1”’ at the clock value plus the sum of those initial times. 
(“PL1” or “PL2”, which “PLProcess” scheduled, changes the value of the target number 
indexed boolean, PLConfirm, to “True’”.) When the “Gun1” is taken off the event queue, 
it checks to see if PLConfirm is “True”, which it will be. Since the PL has confirmed the 
mission, “Gun1” schedules a “Launch”, taking into account the current ranges to the 
target and missile flight times (see “BtryProcess”’ discussion of “GetLaunchTime” 
procedure above). 

The second sub-case occurs if the PL confirmation arrives after the battery 
has finished initial processing of the mission data. (*PLProcess” schedules a “PL1” ora 
“PL2” and a “GunProcess”’ at a time after the time “GunProcess” and “Gun1” finishes 
initial data processing.) In this sub-case, when “Gun1” checks the boolean PLConfirm it 
is “False”. (“PL1” or “PL2” has still not executed to change the boolean value.) 
Therefore, the sub-routine ““GUN1” must wait for the PL confirmation to arrive, turns the 
value of a target number indexed boolean DataProcessed to “True’’, and does not 
schedule a launch. When the “GunProcess” scheduled by “PLProcess’”’ is executed, 
DataProcessed is ““True’’, and a launch is schooled through the procedure 
“GetLaunchTime”. 

The second possible case occurs if the battery does not receive the 
AFATDS data before the PL confirmation arrives. In this case, the PL transfers the data 
to the battery and the battery begins immediate processing of the mission. In the 
simulation, this means that the ““GunProcess” scheduled by “AFATDS” occurs at a time 
after the ““GunProcess” scheduled by “PLProcess”. As the “GunsProcess” scheduled by 
“PLProcess”’ is executed, the value of ProcessingData is “False”. This now becomes a 


PL initiated CFF, and the sub-routine changes the value of ProcessingData and 
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PLInitiated to “True”, collects a total sample of the three lognormal processing times, 
and uses “GetLaunchTime” to schedule a launch. When the “AFATDS” scheduled 
‘““GunProcessed” is eventually executed, the booleans ProcessingData and PLInitiated 


are both ““True’’, so no mission is processed. 


é. “BDEI”’ 


“BDE1” is scheduled by “BDEProcess” and represents the time at which 
the BDE establishes contact with the PL and passes the mission confirmation and firing 
status “Fire” to the PL. The observer and target number indexed boolean variables 
BDEHF and BDEConfirm are changed to “False” and “True”, respectively, indicating 


that BDE wants the mission to be fired. 
Z “BDE2” 


“BDE2” 1s scheduled by ““BDEProcess” and represents the time at which 
the BDE establishes contact with the PL and passes the mission confirmation and firing 
status “Hold Fire” to the PL. The observer and target number indexed boolean variables 
BDEHF and BDEConfirm are both changed to “True”, indicating that BDE does not 


want the mission to be fired. 
g. so ol Ae 5b) 


‘“PLI” 1s scheduled by “PLProcess” and represents the time when the PL 
establishes contact with the battery and passes the mission confirmation and firing status 
to the battery. The observer and target number indexed boolean variables PLConfirm 
and PLHF are both changed to “True”, indicating that the PL has confirmed the mission 


and is passing a “Hold-Fire” to the battery for the mission. 
h. ge ed Sd 


“PL2” is scheduled by “PLProcess” and represents the time when the PL 


establishes contact with the battery and passes the mission confirmation and firing status 


to the battery. The observer and target number indexed boolean variables PLConfirm 
and PLHF are changed to “True” and “False”, respectively, indicating that the PL has 


confirmed the mission and is passing a “Fire” to the battery for the mission. 


I. “GUNI” 


“GUN1” is scheduled by “GunProcess” and represents the time when the 
battery has completed initial processing of the mission data from AFATDS. First, the 
observer and target number indexed boolean variable DataProcessed is changed to 
“True”. The value of the variable PLConfirm is then checked. If PLConfirm is “True” 
and PLHF is ‘False’, “GetLaunchTime” is called to schedule the launch of a missile for 
the target. If PLConfirm is “False’’, nothing further is scheduled, and control returns to 


MASTERTIMER2. 


5. Final Common Subroutines 


a. “Launch” 


The “Launch” sub-routine replicates the “pulling of the trigger” by the 
gunner (State 7 + 4). A U(0,1) draw is made to determine if the missile is lost at launch 
(State 7 +5), misfires and the gunner re-fires (State j + 6), or if there is a successful launch 
(State 7 + 8). The path determined by the Monte Carlo draw determines which sub- 
routine is scheduled, and then returns control to MASTERTIMER. 


b. “Mis Fire” 

This sub-routine replicates the events in State 7 + 6. As discussed in the 
previous chapter, there are four possible transitions from this state. A single U(0,1) 
sample is drawn, which determines the state to which transition is made. If the Monte 
Carlo draw indicates a transition to State 7 + 5 (EFOGM Lost), the number of remaining 


EFOGMSs is decremented and control returns to MASTERTIMER or MASTERTIMER2. 


If the transition indicates another misfire, the ‘“MisFire” sub-routine is scheduled again 
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after a user defined delay and control returns to MASTERTIMER or MASTERTIMER?2. 
A transition to State 7 + 7 (Dud Missile/Select New MissilelReshoot) causes the 
‘“DudRefire”’ sub-routine to be scheduled after a user defined time delay, the number of 
remaining EFOGMs is decremented, and returns control to MASTERTIMER. The last 
option leads to the scheduling of “Flight” (State 7 + 8, Successful Launch) with no delay 


in time, and returns control to MASTERTIMER or MASTERTIMER?. 


c. “DudRefire”’ 


This sub-routine represents State 7 + 7 (DudRefire). It draws a U(0,1) 
random number to determine the succeeding state. Since State j + 7 can transition to the 
same states as discussed in the “Launch” sub-routine, “DudRefire” executes the same 
algorithms, except with different user defined state transition probabilities. It then 
decrements the number of available EFOGMs, and returns control to MASTERTIMER or 
MASTERTIMER2. 


d. “rig? 


This sub-routine considers all possible occurrences during the EFOGM 
flight after a successful launch (State 7+ 8). It calculates the flight time of the missile to 
the target and range of the launcher to the target using the functions TIMEOFFLT and 
RANGE2. (TIMEOFFLT is a procedure that calculates the time the missile will need to 
fly to the target, taking into account the forward velocity of the target. RANGE2 is a 
procedure that calculates the current range from launcher to target based on initial range 
and elapsed time and target velocity.) A U(0,1) draw is first made to determine if the 
missile is detected while in flight. If comparison with a user defined probability of 
detection results in the missile not being detected (State 7 + 10), another U(0,1) is drawn 
to determine if the missile becomes lost while in flight (State 7 + 13). If the missile is not 
lost in flight, transition to State 7 + 14 (EFOGM Target Acquisition/Lock On) is 
simulated, and the sub-routine “Impact” is scheduled to occur at a time equal to the 


current clock value plus the time of flight. Control then returns to MASTERTIMER or 
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MASTERTIMER2. If the missile is lost in flight, the event is considered to have occurred 
at a time value equal to the clock value plus one half of the time of flight value. 

If the missile is detected in flight, another U(0,1) draw is made to 
determine if the enemy destroys the missile in flight (State 7+ 11). If the missile is 
destroyed, the sub-routine ends and control returns tt MASTERTIMER or 
MASTERTIMER? with no further events scheduled. If the Monte Carlo draw indicates 
the missile is not destroyed in flight after enemy detection (State j + 12), a final U(Q,1) 
draw is made to determine if transition is made to State j + 13 orj + 14. Results of these 


transitions were discussed in the previous paragraph. 
é. “Impact” 


The “Impact” sub-routine covers States 7 + 15 andj + 18. A U(Q,1) draw 
determines if the missile impacts on the target (State j + 15) or is a near miss (State j + 
18). If the Monte Carlo draw indicates an impact on the target, an additional U(0,1) draw 
is made to determine if the target was destroyed (State j + 16), or only damaged by the 
resulting explosion (State j + 17). If the second draw indicates the target was only 
damaged, a third U(0,1) draw is compared to user specified probabilities that determine 
the probability of certain damage levels (States 7 + 19 through j + 22), given impact on 
the target. This third draw determines to which damage level state transition is made. 
Once in the damage level state, the damage level to the target is assigned to an array. 
Possible damage levels are “Destroyed”, “Heavy”, “Moderate”, “Light”, and “None” 
(indicating a dud warhead round, or missile impact so far away from the target that no 
collateral damage was caused by the warhead explosion). For simulation purposes, 
damage levels are assumed to be consistent across all enemy vehicle types. (These can be 
readily changed to support user specific analyses.) Once the damage level has been 
changed, the sub-routine “BDA” is scheduled and control returns to MASTERTIMER or 
MASTERTIMER2. 
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If the missile misses the target (State 7 + 18), a U(O,1) draw determines the 
level of damage achieved from a warhead detonation near the target, and determines the 
corresponding transition state. Again, once the damage level annotated, the sub-routine 


“BDA” is scheduled, and control returns to MASTERTIMER or MASTERTIMER?. 


f. “EDA” 


It is assumed that once impact occurs, the FO immediately begins the BDA 
process. The final sub-routine, “_RDA”’, represents the FO’s determination of the level of 
damage caused by the EFOGM missile to a target. The sub-routine searches for a match 
between the target's damage level and the damage levels in the sub-routine. Once a match 
is made, a U(0,1) draw determines if the FO decides to repeat the fire mission on the 
target. The comparison is based upon user defined probabilities of repeating a CFF, given 
the specific damage level to the target. If the U(O,1) draw indicates that a repeat of the 
fire mission is not attempted, the sub-routine ends and contro] returns to the 
MASTERTIMER or MASTERTIMER2. However, if the U(O,1) indicates that a repeat 
CFF is made, “CFF” is scheduled to occur at the current clock time plus a user defined 
amount of time needed by the FO to asses battle damage, and control returns to 


MASTERTIMER or MASTERTIMER?2. 


D. REPLICATION COMPLETION 


The program will continue to schedule events until the battle clock reaches 200 
minutes. At this time, all enemy vehicles moving at the constant rate of 5.55 meters per 
second would be “on top of” the friendly battle position. When this occurs, the 
MASTERTIMER or MASTERTIMER?2 procedure will continue to remove the events 
“Search” and “CFF’, but since either no missiles or FOs are remaining, no further 
processing can occur, and no processing events are added to the Calendar. When the 
time limit 1s reached, program control returns to the main routine that begins to write the 


output data for the replication to a text file. 
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|, DATA COLLECTION 


Within each sub-routine, counters are established to determine the number of 
“visits” to that specific sub-routine. The counter variable VisitTimes[State#,Obsvr,Tet] 
counts the number of times a CFF for a specific FO and target is processed in various 
routines in the model. These routines represent various states in the state space diagram, 
so the number of times a FO acquires a specific target (done in the “Search” sub-routine) 
can be documented. Additionally, time recorders are employed to capture the time each 
event occurs. The variable OutData[ VisitTimes,Obsvr,Tet,State#] records the time at 
which a numbered CFF for a specific FO and target “visits” a specific routine in the 
program. Because of repeat fire missions that may occur, the same FO and target may 
visit the “CFF” procedure more than once. Collection of these data allows comparison of 
the amounts of time required to traverse from state i to state 7. These data are written to a 
text file that can be loaded into a statistical analysis package for use in determining 
distributions and data tendencies. The following chapters discuss the plan for data 


collection and its analysis. 
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IV. ANALYSIS 


A. GENERAL 


The purpose of this chapter is to demonstrate the performance of the simulation 
model within specified parameters, and to demonstrate its sensitivity to any changes in the 
input parameters. The base case for the model’s execution calls for 69 high value targets 
to be engaged by 48 EFOGM. This scenario 1s based upon the simulation depicted in 
HRS 33.5 that was mentioned in previous sections. The mean time and variance values 
used for producing probabilistic time values within the simulation model were taken from 
the Computer Science Corporation (CSC) Technical Report [Ref. 2, Figures 2-4/13] on 
C’ time delays within the EFOGM CFF process. However, during the stress of combat 
and continuous operations, the mean times stated for performance of human tasks by CSC 
in their technical report were considered to be slightly optimistic, based on observations of 
human performance during continuous operations in combat and high intensity training 
environments, such as the National Training Center. Additionally, the processing time for 
the AFATDS computer [Ref. 2, p. 27a] to accept, logically process, and forward the 
numerous CFFs present on the modern battlefield is relatively small, considering the 
abilities of current micro-processors that are able to withstand the durability requirements 
of DOD hardware and the environmental extremes that may be met in combat operation 
areas. Accordingly, mean times for human tasks (1.e., CFF processing, intercept 
prediction plotting, issuing firing orders) for both the digitized and non-digitized cases, 
and the mean processing time for the AFATDS computer in the digitized case, were 
increased by 25% to reflect the perceived time underestimation. Simulation runs executed 
under these new mean time values constituted the base case + 25% replications. Values 
for the variance for each adjusted mean were assumed to have remained the same. 


Replications were made for both the digitized and non-digitized scenarios using the two 
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input time parameters. The following sections discuss the base case and the base case + 


25% data collection and replication outcomes. 


B. DATA COLLECTION 


For each replication run, time events (measured in minutes) were recorded for use 
in data analysis. Time events in both the digitized and non-digitized cases were recorded 
for a specific target and observer. Similar events for the digitized and non-digitized 


replications include: 


e Time of detection (T1). 

e Time the battery completes CFF processing (T8). 

e Time a missile was lost at launch (LAL), failed to fire or is a DUD round (T9). 
e Time of a successful launch (T10). 

e Time of impact (T11). 


Events specific to the non-digitized replications include: 


e Time the CFF arrived at the brigade HQs (N2). 

e Time the CFF arrived at the FSO (N3).. 

e Time the FSO began processing the CFF (N4). 

e Time the CFF arrived at the LNO (NS). 

e Time the LNO began processing the CFF (N6). 

e Time the battery began processing the CFF (N7). 


Events specific to digitized replications include: 


e Time the CFF arrived at the AFATDS for processing (D2). 

e Time the CFF arrived at the BDE cell (D3) after being sent from AFATDS. 
e Time the CFF arrived at the PL cell (D4) after being sent from AFATDS. 
e Time the CFF arrived at the Battery (D5) after being sent from AFATDS. 

e Time the PL receives mission confirmation from the BDE cell (D6). 

e Time the battery receives mission confirmation from the PL (D7). 


To compare/contrast the differences between the digitized and non-digitized systems, the 
total CFF processing time was needed. The total CFF processing time is simply the 


difference between the time the target was detected (T1) and the time the battery finishes 
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processing the CFF (T8). The output matrices of event times organized by observer and 
target were then converted into spreadsheet format. (See Appendix N for an example of 
an output data file.) Simple column operations determined the time difference between T1 
and T8 for each CFF that resulted in a missile launch. Twenty-five simulation replications 
produced 1200 values for total CFF processing time. Additionally, for the non-digitized 
case, the time a CFF spent “waiting” in the FSO or LNO queue was determined by finding 
the difference between event times N4 and N3 (FSO queue time) and event times N6 and 
N5 (LNO queue time). These columns of numbers were then used to create descriptive 
Statistics and estimate probability distribution fittings for the processing times produced by 


replications of the simulation. An analysis of these times will now be discussed 1n detail. 
C. NON-DIGITIZED REPLICATIONS 


1. Base Case 


Using 25 non-digitized replications of the simulation model, the descriptive 


statistics contained in Table 2 were obtained. 


State 8 Time - State 1 Time = | Non-Digitized Data 
Processing Time 









ie on 2 Saal 9,328 
Standard Error 0.099 
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Table 2. Descriptive Statistics, Non-Digitized Replications, Base Case 
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The histogram shown in Figure 4 shows the shape of the data for the processing times of 
the non-digitized base case replications. 


NON-—DIGITIZED DATA HISTOGRAM 
25 REPLICATIONS = 1200 SAMPLES 


100 150 200 


NO OF SAMPLES 
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Figure 4. Histogram for Non-Digitized Data 


Using the program BestFit©, the 1200 data points were analyzed and specific distributions 
ranked according to the chi-square, Kolmogorov-Smirnov and Anderson-Darling 
goodness-of-fit tests. This allowed the elimination of unlikely distributions from 
consideration. Choosing only the top three distributions for extensive analysis, the 
statistical package AGSS was used to specifically determine how well the data followed 
the chosen probability distributions recommended by BestFit©. The exploratory data 
analysis conducted by AGSS identified the data came nearest to matching the lognormal (U 
= 2.168, 6 = 0.3563) or Gamma (a = 7.916, B = 1.79) distributions. Graphical displays 


showing the density function, cumulative distribution function, cumulative hazard 
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function, the probability plot and the analysis for the lognormal distribution fit are shown 
in Appendix J. While neither the Kolmogorov-Smirnov, Cramér-Von Mises, nor the 
Anderson-Darling goodness-of-fit tests demonstrate an exceptional fit, the graphical 
representations indicate the two distributions provide a relatively good match over the 


range of the data. 


2. Base Case + 25% 


Executing the simulation model using the non-digitized replication runs, and with 
time values adjusted as discussed in Section A, the data in Table 3 were obtained. 
Examining only the difference in the mean processing time for the base case and base case 
+ 25% shows an increase in total processing time of over 8 minutes. The increase in the 
base case + 25% mean processing time 1s not a straight 25% increase over the base case 


mean processing time. The causes of this increase will be discussed later in the chapter. 


State 8 Time - State 1 Time =| Non-Digitized Data 
Processing Time 
ti tei Se 


-0.938 


5.844 


Maximum 35.733 
sumo 20591.111 
1199 
35.733 


Smallest(1 5.844 


Confidence Level(95.000% 0.393 





Table 3. Descriptive Statistics, Non-Digitized Replications, Base Case + 25% 


The histogram shown in Figure 5 shows the shape of the data for the processing times of 


the non-digitized, base case + 25% replications. 
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Figure 5. Histogram for Non-Digitized Data, Base Case + 25% 


Using BestFit©, the empirical data was examined for a possible fit to probabilistic 
distributions. The two distributions chosen by BestFit© were examined by AGSS which 
showed the data followed the Weibull (a = 2.708, B = 19.292) and Beta (a = 2.729, o 
= 2.924) distributions. AGSS showed no significance within the standard goodness-of-fit 
tests, and graphical analysis showed the empirical distribution of the data does not closely 


follow either distribution. (See Appendix K). 
D. DIGITIZED REPLICATIONS 


1. Base Case 


Again, 25 simulation replications were executed to generate 1200 data points for 
the digitized base case version of the model. Shown in Table 4 are descriptive statistics 


for the processing time for those replications. 
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State 8 Time - State 1 Time = Digitized Data 
Processing Time 


Mean si sis 1.456 
0.004 
Median sid 1.447 

ode 1.447 
0.018 
2.169 
0.829 
Range ——is—s—s—s 1.098 
Minimums 
Maximum sid 2.218 
Sum ——sSCisS 


um 1747.67 


1200 


argest(1 2.218 


mallest(1 
Confidence Level(95.000% 0.008 


Table 4. Descriptive Statistics, Digitized Replications, Base Case 
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Shown in Figure 6 is a histogram showing the general shape of the data for the digitized 


base case CFF processing time. 
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Figure 6. Histogram for Digitized Data, Base Case Replications 
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Using BestFit© and AGSS the empirical data was shown to fit the lognormal (u = 1.456, 
Oo = 0.0889) distribution at a very reasonable significance level for the Kolmogorov- 
Smimov, Cramér-Von Mises and Anderson Darling goodness-of-fit tests. The graphical 


depiction and statistical analysis of data fit are shown in Appendix L. 


2. Base Case + 25% 


Similarly, 25 simulation replications were executed for the digitized, base case + 


25% version of the model. Descriptive statistics for the processing time are contained in 


Tables: 
State 8 Time - State 1 Time = Digitized Data 
Processing Time 


Mean gl ee 
Median 1.743 
Mode 8D 


Sample Variance 0.019 
Kurtosis 1.630 
kewness 0.504 
1.284 
1.385 
2.669 

2097.625 


Table 5. Descriptive Statistics, Digitized Replications, Base Case + 25% 












| 


The empirical distribution of the 1200 processing time data values from these replications 


is Shown in the histogram in Figure 7. 
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Figure 7. Histogram for Digitized Data, Base Case + 25% 

Similar to the base case for the digitized replications, BestFit© showed a strong 
correlation with the lognormal! distribution. AGSS confirmed this finding at a high 
significance level for all goodness-of-fit tests, and showed the data to fit the Lognormal 
(uu = 1.748, o = 0.138) Distribution. Graphical and statistical descriptions of the 
distribution fitting are shown in Appendix M. Contrary to the increase in the mean 
processing time between the Base Case and Base Case + 25% cases in the Non-Digitized 
replications, the mean processing time in the Base + 25% case of the digitized replications 
is almost exactly 25% higher than the base case, indicating the increase in mean processing 
time to be purely a function of the increased mean values of the input parameters. 

Another concern would be to check if the digitized system remained constant over 
the range of target availability. To test this premise, both digitized scenarios were 


executed with the standard 69 targets, and then with 35, or approximately half the number 
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of original targets available. Unsurprisingly, the digitized processing times remained 


constant within the two cases over the range of targets available. 
E. COMPARISON/ANALYSIS 


As shown in the previous sections, there is a marked difference in the performance 
of the digitized and non-digitized ec system for the EFOGM. However, this performance 
is based on a fixed number of targets, moving at a constant march rate, and at 
predetermined intervals, arriving into the target area. To determine the sensitivity of the 
model to changes in the number of targets available to the EFOGM C° system, a series of 
replications were executed with decreasing numbers of targets available for prosecution. 
The decreasing numbers were accomplished by “thinning” the target array by five targets 
at atime. Targets were not all taken from the middle or rear of the target array, but 
evenly removed from the array of high value targets approaching the battle area. This has 
an effect similar to adjusting the arrival rate of the targets. The ultimate goal is to 
determine at what point, in the number of targets approaching the battlefield, does the 
FSO and LNO queue become overloaded enough to significantly reduce the efficiency of 
the processing system. The premise behind this goal is that the queue times in the non- 
digitized system heavily influence the total processing time for the system. However, if 
the number of presented targets is such that one target can be acquired and processed for 
prosecution before the next target becomes available, then the FSO and LNO queue are 
not a factor in the total time needed to prosecute the target. Through the process of 
“divide and conquer’, the number of targets available was halved, then halved again, and 
so on until the mean times in the FSO and LNO queues became approximately zero. At 
this point, processing time can only be attributed to the amount of time required to 
process targets, and not a function of time waiting to be processed. Shown in Figure 8 is 
a composite graph depicting the increase in the mean base case processing time for non- 
digitized replications versus the mean base case processing time for digitized replications. 


The distribution of results for the non-digitized case are further represented by box plots. 
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DIGITIZED vs NON-DIGITIZED PROCESS TIMES 
BASE CASE, 10 REPLICATION RUNS PER SAMPLE 
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Figure 8. Composite Graph, Digitized/Non-Digitized Target Processing Times, Base Case 


From about [5 to 35 targets, the increase in the mean and variance of the 
processing time for non-digitized replications remains relatively small, and the delta 
between digitized and non-digitized CFF process time remains almost constant. However, 
at 40 targets, the mean and variance for the corresponding target numbers begin to show a 
marked increase, while the processing time for the digitized replications remains constant. 
Further analysis indicated that for more than 40 targets in the target array, new targets 
arrive in the area faster than the non-digitized system can process previously arrived 
targets. Therefore, new targets are placed on either, or both the FSO and LNO queues 


before processing. Shown in Figure 9 is a graph depicting the overall non-digitized base 
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case CFF processing mean times for increasing target numbers, and the corresponding 


mean time the CFF spent in the FSO and LNO queue. 
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Figure 9. Comparison of Non-Digitized CFF Processing Queue Times, Base Case 


Figure 9 shows that the difference between the LNO processing time and the Mean Total 
CFF Processing time remains basically constant over the range of targets. It is, in fact, the 
FSO and LNO queues which cause the value of Mean Total CFF Processing time to 
increase. 

In the non-digitized base case + 25% replications, the mean CFF processing times 
are significantly larger than in the non-digitized base case. Shown in Figure 10 is a 


composite graph of the total mean processing time of the non-digitized, base case + 25% 
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replications, with comparison to the digitized case. The distribution of times in the non- 
digitized case are further represented by box plots. The absence of box plots for the 
digitized cases indicate the surplus capability of the digitized system. 

DIGITIZED vs NON-DIGITIZED PROCESS TIMES 

BASE CASE + 25, 10 REPLICATION RUNS PER SAMPLE 
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Figure 10. Composite Graph ,Digitized/Non-Digitized Target Processing Times, Base Case + 25% 


In this case, with approximately 30 targets in the target array, the CFF processing time 
begins to increase dramatically for both mean and variability, showing the increasing tails 
of the box plots and the appearance of more outliers. As in the non-digitized base case, 
as the numbers of targets available for prosecution increase, the mean time a CFF spends 
in either, or both, the FSO and LNO queues dramatically increases, thereby increasing the 


total mean processing time. This result is shown in the graph in Figure 11. 
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Figure 11. Comparison of Non-Digitized CFF Processing Queue Times, Base Case + 25 % 


Another question is, “If the number of targets that are available to be prosecuted 
remains constant, what would occur in a case where there 1s a greater supply of missiles 
than targets?” For this excursion from the base case, the number of targets was returned 
to its original value of 69 and the number of missiles was increased to 100. This change 
allowed more re-engagement of targets that were not destroyed by the first missile 
processed for that target. In the base case, the mean time a CFF spent in the FSO queue 
increased almost five minutes (1.6 to 6.5 minutes); time in the LNO queue increased over 
three minutes (3.9 to 7.3 minutes); mean total processing time increased from 9.9 to 15.1 
minutes. For the base case + 25% the increases were even larger: 4.1 to 13.4 minutes in 


the FSO queue, 7.2 to 11.9 in the LNO queue, and 15.8 to 21.9 in total mean CFF 
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processing time. While the increase in the number of available missiles means more 
missions can be fired, processing the CFFs for those missions swamps an already 
overburdened non-digitized system. Hence, the drastic increase in the CFF processing 
times for both the base case and base case + 25% scenarios. However, for both the base 
case and base case + 25% digitized scenarios, the processing time remained constant (1.45 
minutes in the base case, 1.75 minutes in the base case + 25%), clearly indicating that the 
digitized system has not yet reached capacity. A matrix of all runs with mean times and 


standard deviations is located in Appendix N. 


in 
Cn 








V. CONCLUSIONS 


ee SUMMARY 


While the simulation model developed here has some limitations, there are some 


definite benefits to such discrete time simulations: 


e The overhead cost is inexpensive. Pascal is a low-priced software program, and 
not difficult to learn. [It runs on standard DOS based computers currently present 
in an overwhelming majority of Army analysis agencies. The code is easily 
understood, even by those who have little experience with computer programming. 


° The discrete time simulation approach provides a quick, reasonable answer to an 
existing problem. With minor adjustments to the programming code, the model 
easily becomes highly robust (see discussion below). 


° The outcome from these initial simulation runs provides a firm basis for future and 
follow-on studies concerning the nature of the problem and possible changes to the 
system before production. 


Examining the data presented in the previous section, several conclusions can be 


reached: 


e The digitized system is far superior to the non-digitized system regardless of the 
number of targets. 


e The non-digitized system’s processing remains constant when target movement 
rates are slow (approximately 10 meters/second) or when a limited number of 
targets are present (approximately 40). 


° Increasing the supply of missiles has a significant effect on non-digitized queue 
times. With fewer missiles than targets, after the first round impacts on a target 
and BDA is performed, there generally were no more missiles remaining to re- 
shoot the target if it was not completely destroyed. However, with more missiles, 
a CFF can be reinitiated on the target, thereby increasing the number of CFFs on 
an already overburdened system. 


° The time a CFF spends in the FSO/LNO queues has the greatest impact on total 
CFF processing time. With more targets to process, increases in individual 


oe, 


processing times, or more missiles than targets, the queues become bottlenecks 
that continue to accumulate CFFs. 


e With the AFATDS being capable of processing multiple CFFs, and simultaneous 
CFF processing occurring at the brigade, platoon leader and battery levels, the 
total processing time is purely a function of how fast data can be manipulated and 
disseminated. 


° The stress of continuous combat operations may have a significant effect on human 
processing abilities as seen in the non-digitized base case + 25% scenario. This 
stress could greatly impact total processing times for CFFs. 


These conclusions are based upon the results of the simulation model which in 
itself has not been validated nor verified. Further testing of the model and its parameters 


would need to be made to come to more verifiable results. 


B. LIMITATIONS AND ENHANCEMENTS 


There are several limitations to the discrete ttme model. While they can be 


overcome they did cause some problems during the development of this thesis. 


e Limitation of the Pascal language make programming difficult at times. With a 
limited allowable data segment (65,520 bytes), this restricts the number of 
variables that can be declared within a procedure. With the size of the data arrays, 
the data segment limitation can easily be reached. This causes the programmer to 
parse out global variables to several sub-unit programs, and can make the code 
difficult to follow at times. 


e The programming language is not object-oriented in nature. While object oriented 
programming is possible in Pascal, it is not as powerful a language as C™ or 
MODSIM would be for developing this type of program. However, C and 
MODSIM can be difficult to understand, and MODSIM, while extremely 
powerful, is usually run on UNIX based systems that have substantially higher 
costs than normal DOS based hardware applications. PC versions of MODSIM 
exist, but are very expensive. 


e Pascal can be very difficult to debug and remove run-time errors. The Pascal run- 


time library provides error messages, but these messages can be deceiving and 
difficult to trace, especially over numerous sub-programs and routines. 
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° No queue is depicted within the digitized C’ system. While the new AFATDS 
system may be powerful, it is doubtful it can manage the volume of CFFs that 
could theoretically be placed on it. Additionally, there is no system of checks and 
balances depicted in the CSC diagrams. In a force that is seriously concerned with 
the possibilities of fratricide on the battlefield, not having some command and 
contro] headquarters wired into the decision loop for executing missile firing 
criteria is an error. The simulation model seeks to correct this by allowing no CFF 
to be executed without the tacit approval echoing down the chain of command to 
the battery. In the “do not fire by exception” policy, the battery could theoretically 
receive a fire mission from the AFATDS and execute that mission without the 
brigade cell or platoon leader ever receiving the firing data from the AFATDS. 
Should some command and control issue arise at brigade that would restrict 
missile flights (1.e., attack helicopters engaging deep targets), there are no 
provisions shown in the CSC report calling for explicit brigade approval for missile 
flights. 


° Only one “super-launcher”’ is depicted in the simulation model. In reality, there 
would be several launchers situated in numerous locations over the friendly battle 
position to enhance survivability. This would require different ranges from each 
launcher to the target array, and cause the possibility of communications failures to 
occur while passing fire commands from the battery to gun level. Also, while it 
would seem that counter-battery fires would be difficult on these launchers, no 
launcher attrition was depicted, nor were the possibilities of “shoot-and-move” 
style tactics taken into consideration. These tactics would obviously render 
launchers inactive for the period of time they are moving, and unable to shoot. 
(This also raises the possibility of hip-shoots.) 


° Finally, enemy targets were depicted as approaching the battlefield arrayed in a 
single file march column formation with doctrinal distances between each vehicle 
and maneuver element within the attacking force. In reality, the enemy will most 
likely attack on a broad front across numerous axes of advance, which only 
increases the command and contro] burden on the EFOGM system. 


C. VARIATIONS AND MODEL CAPABILITIES 


While there are some obvious limitations to the model, there are some inherent 
capabilities and variations of model runs which make this model highly robust and 
versatile. This thesis concentrated on only one aspect of the EFOGM system, the 
processing time for the two possible C° systems. However, several features of the model 


are worth mentioning: 


Spy 


e Missile characteristics can be adjusted and studied. The model allows for missile 
detections by enemy forces and evasion. While the probability of occurrence of 
these events in based on the author’s professional judgment, “exact” probabilities 
based on other studies can easily be input to give estimates of missile successes 
while considering enemy air defense threats. Additionally, the probability of target 
damage based on a missile hit, the probability of a missile hit, and several other 
missile characteristics were all based on professional judgments of what they might 
be. Actual data from AMSAA or TEXCOM would only increase the validity of 
the model. 


e Detections by sensors, while realistic, are based on relatively old detection theory 
formulas. New detection algorithms exist which would allow the detection portion 
of the simulation to be much more realistic. For example, the NVEOL (Night 
Vision Electro-Optical Laboratory) detection algorithm simulates the capabilities 
of image intensification and thermal sights which are not taken into account in this 
model. Target detections based on a more sophisticated detection algorithm could 
easily be fed to the simulation model, thereby increasing the realism. 


° As the C° system for the EFOGM is further refined through testing new and more 
highly advanced technology, the simulation model provides a framework to test 
new capabilities for different outcomes. 


This thesis demonstrates the potential of relatively simple, discrete-time simulation 
models to assist in performance estimation of undeveloped combat systems on the 
modern-day battlefield. Building on the principles of the stochastic state space, it provides 
analytic agencies, specifically, the U.S. Army Infantry School, the study sponsor, the 
ability to incorporate future technologies into existing combat models, or be used as a 
stand-alone model for data analysis. In light of today’s quickly developing technology 
base, and the unpredictable nature of future conflicts, simulation modeling can only 
improve the ability of Army planners to prepare for future missions, and strengthen their 


commitment to provide the Army with the best equipment available. 
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APPENDIX A. DIAGRAM EXTENDED CLOSE BATTLEFIELD 


| 83.7 Kilometers | 


Max Range of 
FO/Sensor -10 km 


2km - Min 
Range of 
EFGOM 


EFOGM 
Battery 





Figure 12. Diagram of Simulated ECB 
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APPENDIX B. STATE SPACE TRANSITION PROBABILITIES 


The transition probabilities are based on the author’s professional judgment. Refer 


to Figure 1 for graphical depiction of these state space descriptions. 


| ft | FOinitialization | So 2or5 | .015,0.985_ 
a 


5 FO Undetected and — 
not Destroyed 


FO Searches 7 through j or 1 Should a detection be made, the 
transition will be made to the type 
of detected system in states 7 
through j with probability = 1.0. 
However, if a detection is not made 
in the described time period, 
transition is back to state 1 with 
probability = 1.0. 


_ ON ee ee a re 
MP DetectADA asset | eC 
Ma) |) etectjthiassets |e OT 


Jtl FO transmits a CFF jt+2o0rj+3 ON ie 0. 99 
to the BDE or 
AFATDS 


CFF not received by jt+t2orj+3or] The probability the FO makes 
BDE or AFATDS contact with BDE/AFATDS on his 
next try is 0.95. Chance he fails to 
make contact again is 0.05. If the 
number of attempts cause the time 
standard to be exceeded, the FO 
quits, and returns to State 1 with 
probability = 1.0 

































j+3 CFF is received by 
BDE/AFATDS and 
processed by all 


subordinate elements 


| j+4 | EFOGM launch 0.01, 0.02, 0.97 
jt+5 EFOGM lost at None N/A 
launch (LAL) 


| j+6 | EFOGMMisfires | (+5, )+6,1+7 orj +8 0.01, 0.05, 0.01,0.93 
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Transition Probability (Pij) 


Jt+7 DUD missile, select j+5,j+6j+8 0.01, 0.02, 0.97 
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APPENDIX C. RANDOM NUMBER GENERATOR 


This appendix contains the PASCAL code used within the simulation model to 
generate random U(0,1), normal and lognormal numbers. The U(0,1) random number 
generator is documented in the Law and Kelton text [Ref. 7, pp. 451-454]. The other 


random number generators were developed for this study in accordance with the text. 
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{ AUTHOR: CPT David S. Pound 


ASSIGNMENT: Thesis 
WRITTEN: 9 April 1994 
REFERENCE: Simulation Modeling and Analysis, 2d Ed.,Law and Kelton, pps.451-54. 
OBJECTIVE: A random number generator for thesis simulation program with prime 
modulus multiplicative linear congruential generator: 
Z[i)] = ( 6380360016 * Zli=1) (MOD 2147483647), 
based on Marse and Roberts'portable FORTRAN random-number generator UNIRAN. 
Multiple (100) streams are supported, with seeds spaced 100,000 apart. 
Throughout, input argument Stream must be an Integer giving the desired stream 
number. The initialization procedure RANDDF described below must be invoked 
before using the generator, in order to set the seeds for the 100 predefined streams. 
Usage: (Four procedures) 
1. Before using the generator, it is required to initialize the routines by 
executing 
Randdf ; 
This sets the initial seed values for all 100 streams in the array Zrng. 
2. To obtain the next U[0,1] random number from stream STREAM, execute 
U := Rand(Stream) 
The Real variable U will contain the next random number. 
3. To set the seed for stream STREAM to a desired value ZSET, execute 
Randst (Zset, Stream) 
where Zset must be an Integer constant or variable set to the desired seed 
a number between 1 and 2147483646 (inclusive). Seeds for all 100 streams 
are given in the code, and must be initialized by invoking Randdf. 
4. To get the current (most recently used) integer in the sequence being 
executed for stream STREAM into the Integer variable Zget, execute 
Zget = Randgt (Stream) 


Unit RandNumGen; {Ref: Simulation Modeling and Analysis, Law & Kelton, pp 451-54} 
INTERFACE 


VAR 
Zrng : ARRAY [{1..100] of LongInt; 
Procedure RANDDF; 


Function RAND(Stream : LongInt) : Real; 

Procedure Randst(Zset, Stream : LongInt); 

Function Randgt(Stream : LongInt) : LongInt; 

Function Seed: Word; 

Function Normal (Mu,Sigma :Real) : Real; 

Function LogNormal (Mu, Sigma :Real) : Real; 

IMPLEMENTATION 

uses CRT, DOS; 

{------------------------ Set the seeds for all 100 Streams-------------------- } 

Procedure RANDDF; 

Begin 

Zrngl 1) :=19732Z7/2912: Zrng(\ 2): = 281629770 Zrngies) = scuoleZ U. 
Zrngl 4] ;=1280689631:>- Zingl 5)):=2096730329;>" Zrngl ol =1o3se7 culo, 
Zrng{ 7]:= 913566091; Zrngi{ 8):= 246780520; Zrng[ 9] -=1363774876- 
Zrng(10):= 604901985; Zrng[11)]:=1511192140; Zrngf12))-=19335576050- 
2rngi[l3a):— 824064364.  Zrng([14):=— 150493254. Zrng (ls) 2 24270 soe 
Zrng(16)== °752534171;, Zrng{17) ;=1964472944- Zrng ils) =12027299975- 
Zrngll9) <= 233287322; 2rng (20) s=19 11 Zis000} 2rng (Zi = 9725357 oss. 
ZrngizZ):= 42403492145>" Zrngi 23) <= 995232223 2rInGg | 24 a=) 0320S oo 
Zrng(25]:= 762430696; Zrng[i26)]:=1922803170- Zrng iz ajc ole e2s 
Zrng(28):= 76271663; Zrng[29]:= 413682397> Zrngid0)-— 97 264cccu4- 
Zrng(31):= 336157058; Zrng[(32)2=1432650381>. 2Zrnclse) = 046304 
2rng(34)]:= 595778810; Zreng[35]:= 377722390; Zrnglael-—104e5 244—>- 
Zrng(37J:= 68911991; Zrng(38) :=2088367019; Zrng(39]:= 748545416; i 
Z2rng{40):= 622401386; Zrng[(41] :=2122378630; Zrng[4Z)<= 640690903; 
Z2rng(43] :=1774806513; Zrng[44] :=2132545692- Zrngl45)-=2079249570. 
2rng[46):= 78130110; Zrngl47)-= 852776735; Zrng(43l-—=llevecr 272 
Z2rng[{49) :=1351423507; Zrng(50)] :=1645973084; Zrng(51] :=1997049139; 
2rng([S2]:= 922510944; Zrng(53]-=20455128670- Zrng[(S4)e— vel socs7 7 
Zrng [55] = 243649545; Zrng[56] s=100481e772>" Z2rngl571-—8 7 (so ccUG 
Zrng([58):= 403188473; ZrngiS59]+= 372279677 -) 2eng oO b.—- "oO rows 4e5- 
Zrng[61]:= 498067494; Zrng{62]2=2087759556-" 2rnGgl os —.4951> 7 o> 
z2rng{64]:= 597104727; Zrng[65] :-=1530940798- Zrngl66]- =1361 4496274. 
Zrng{67]:= 536444882; Zrng/[68] -=16631S3658= 2Zrng (Gol = 555047 25- 
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Zeng (/0):=> 67784357; Zrng[71)+=1432404475; Z2rng[72]2:= 619691088; 
Zeno sie = 229025595; Zrng (74)2= 8303802310; >z2rng (75) s= 176192644; 
Zeng 7ol-=11167980070;> Z2ing(77)-=. 277854671. \2rng [78] 27=13665980550; 


Beng (79) :=1142483975; 2Zrng[80)]-=2026948561. > Zrng [81] :=1053920743; 
arng | 62) +=" 786262391." Zrng|63)]-=1792203830;-2rng [34 )+=1494667770; 
Zrng (85) :=1923011392; Zrng[86] :=1433700034; Zrng[87] :=1244184613; 
Gro sa)2=1147297105-> Zeng (So) = 539712780.) 2rng|[90]2=1545929719; 


Zrng[91)]:= 190641742; Zrng[(92] :=1645390429; Zrng[93]:= 264907697; 
Zeng o4) = 62Use9255> Zeng losis =1502071852- Zrngl96}2= S27 711160; 
Zrng(97]:= 364849192; Zrng[98] :=2049576050; Zrng[(99]:= 638580085; 
Zrmeg(l00}:= 547070247; 
End; {Randdf} 
aa SS Generate the next Random Number------------------------- } 
Function RAND(Stream : LongInt) : Real; 
CONST 
B2e15 = 32768; 
B2E16 = 65536; 
Modlus = 2147483647; 
Nal e1 = 24112; 
Mult2 = 26143: 
VAR 
Feeko> Hi3sl, LowlS, Lowprd, Ovflow, Zi: Longint; 
BEGIN 


{Generate the next random number} 

Zie = Zrng [Stream]; 

Pao = 71 DIV B2ZE16; 

Beword:= (Z1 - HilS*B2E16) * Multi; 

iewl5:= Lowprd DIV B2E16; 

Beets = Hil5S*Multl + Lowl5: 

Owmetow:= Higl DIV B2E15; 

Zi:= (( (Lowprd-Low15*B2E16)-Modlus) + (Hi31-Ovflow*B2E15)*B2E16) + Ovflow; 
imeaZt <0 THEN © Zisv=-Zi + Modlus: 

hiss = 271 DIV B2E16; 

Mewordad:= (Zi - Hil5*B2E16) *Mult2; 

Low15:= Lowprd DIV B2E16; 

Pasit:= Hil5*Mult2 + Lowl15; 

@yvtlow:= Hi3l DIV B2E15; 

Zi:= (((Lowprd-Low15*B2E16)-Modlus) + (Hi31-Ovflow*B2E15) *B2E16) + Ovflow; 
If Zi < O THEN Zi:= Zi + Modilus; 

Zrng(Stream] := Zi; 

Rand: = 9 (2* (21 DIV 256) + 1)/16777216.90; 

End; {Rand} 


{---------------- Set the current Zrng for stream STREAM to Zset--------------- } 
Procedure Randst(Zset, Stream :LongInt); 
Begin 
Zrng(Stream]:= Zset; 
End: 
ae = Return Che current Zrng fOr stream STREAM----=---=----—---=-==- } 
Function Randgt(Stream : LongInt) : LongInt; 
Begin 
Randgt:= Zrng[Stream] ; 
End; 
{oor ere rrr rr er rere SEED GENERATOR FOR RNG-------------------------- 7 


Function Seed: Word; 
Var Hour,Minute,Second,Sec100: Word; 
Begin 
GetTime (Hour,Minute,Second,Sec100); 
SEED:= Second; 


end; 
{---------------- RETURN A NORMAL (MU,SIGMA) RANDOM VARIABLE-------------------- } 
Function Normal(Mu,Sigma: Real) : Real; 
Var U1,U2,V1,V2,W,X,Y :Real; 
Done : Boolean; 
Begin 
DONE:= FALSE; 
REPEAT 
U1:= RAND(SEED) ; 
U2:= RAND(SEED); 
Vie= 2°Ul-1; 
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V2:= 2*U2-1; 
W:= SOR(V1) + SOR(V2); 
If Ww < 1.0 then begin 
Y:= SORT((=2.0*1n(W))7W) ; 
Done:= True; 
end; {if} 

UNTIL DONE; 

» Ca Bk ae 

Normal:= Mu + Sigma*X; 
End; {Normal} 
{----- eee een RETURN A LOGNORMAL (MU, SIGMA)” RANDOM VARIABLE -———— ee 
Function LogNormal (Mu,Sigma: Real) : Real; 
Var Ul,U2, Vl, V27W. x Y,2,MuzZ,Sigmasar = Real 

Done : Boolean; 

Begin 

Done:=False; 

Mu2:= Ln( SQOR(MU)/ SOQRT( SQOR(SIGMA) + SOR(MU) ) le 

SigmaSqr:= Ln( ( SOR(SIGMA) + SQR(MU) ) / SOR(MU) ) 

Z:= NORMAL(MU2, SORT(SIGMASOR) ); 

LogNormal:= Exp(Z); 
End; {LogNormal } 
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APPENDIX D. VERIFICATION OF RANDOM NUMBER GENERATOR 


The graphs and statistical analysis shown in this appendix were used for 
confirmation of the normal and lognormal random number generators used in the 


simulation program, and were created by AGSS. 


1000 GENERATED NORMAL(0.20, 0.02778) VALUES 


NORMAL DENSITY FUNCTION, N=1000 NORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=1000 
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Figure 13. Normal Graphical Analysis of Randomly Generated Numbers 
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ANALYSIS OF NORMAL DISTRIBUTION FIT 


DATA : NORM 

SELECTION : ALL 

X AXIS LABEL: NORM 

SAMPLE SIZE : 1000 

CENSORING : NONE 
FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKEL|HOOO 
CONF METHOO : EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER — UPPER W SIGMA 
WU 0.24971 0.24801 0.25141 7.5019E-7 0.000E0 
SIGMA 0.02739 0.026253 0.02866 0.0000E0 3.751E-7 
LOG LIKELIHOOD FUNCTION AT MLE = 2178.7 


SAMPLF = FITTED GOODNESS OF FIT TESTS 
MEAN =: «(0.24971 = (0.24971 CHI-SQUARE : 9.3325 
STD DEV: 0.027403 0.02739 * DEG FREED: 7 
SKEWNESS: -0.0097735 0 SIGNIF =: 0.61946 
KURTOSIS: 3.0932 3 KOLM-SWIRN : 0.021278 


+ BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF =: 0.7998 
CRAMER-V M : 0.077831 


PERCFNTILES SAMPLE = FITTED SIGNIF oo: > 15 
9: 0.2058 ~=— 0.20460 ANDER-DARL : 0.4376 
10: 0.21345 = 0.2146 SNIP ats 
29: 0.2324 = 0.23124 
30: 0.2486 = 0.24971 KS, AD, AND CV SIGNIF. LEVELS NOT 
7: 0.268 0.26818 EXACT WITH ESTIMATED PARAMETERS. 
90: 0.28515 0.28482 
95; 0.2949 0.2947 NOTE: A SMALL SIGNIFICANCE LEVEL 


(EG. Ps.01) INDICATES LACK OF FIT 
CHI-SQUARE GOODNESS OF FIT TABLE 


LOWER UPPER 08S EXP = OE_-—((O+E) 2) 
“INF. 0.184 11 8.7232 2.2768 0.59426 
0.1846 0.20138 26 30.104 4.1044 0.55958 
0.20138 0.21816 90 85.884 4.1158 0.19724 
0.21816 0.23495 160 170.22 10.216 0.61318 
0.23495 0.25173 254 234.43 19.568 1.6333 
0.25173 0.26851 216 224.4 8.3992 0.31438 
0.26851 0.28529 145 149.28 4.2813 0.12279 
0.28529 0.30207 65 69.007 4.0071 0.23268 
0.30207 0.31886 27 22.159 4.8412 1.0577 
0.31886 +INF. 6 5.7934 0.20663 0.073696 
TOTAL 1000 1000 5.3325 


Figure 14. Normal Statistical Analysis of Randomly Generated Numbers 


70 


1000 GENERATED LOGNORMAL(0.41667,0.0611) VARIATES 
LOGNORMAL DENSITY FUNCTION, 1V=1000 LOGNORUAL CUMULATIVE DISTRIBUTION FUNCTION, v=1000 
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Figure 15. Lognormal Graphical Analysis of Randomly Generated Numbers 
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ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA : LO 

SELECTION =: ALL 

X AXIS LABEL: LOG 

SAMPLE SIZE: 1000 

CENSORING : NONE 
FREQUENCIES ; 1 

EST. METHOD ; MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 


PARAMETER ESTIMATE LOWER UPPER WJ S1GMA 
MU —).8883 -).89729 +).87932 0.000020995 0 
S1GHA 0.1449 0.13888 0.15162 0 0.000010497 
LOG LIKELIHOOD FUNCTION AT MLE = 1401.1 
SAMPLE = FITTED GOODNESS OF FIT TESTS 
MEAN =: (0.41966 0.4157 CHI-SQUARE : 7.9602 
STD DEV: 0.059929 0.06055 DEG FREED: 8 
SKEWNESS: 0.2737 0.44007 SIGNIF =: 0.43737 
KURTOSIS: 2.9854 5.3463 KOLH-SHIRN : 0.020521 
+ BASEO ON MIDPOINTS OF FINITE INTERVALS SIGNIF =: 0.79367 
CRAMER-V M: 0.048731 
PERCENTILES SAMPLE —- FITTED SINIF =: > .19 
5: 0.32665 (0.32411 ANDER-DARL : 0.32735 
10: 0.3403 0.34163 SIGNIE : aela 
25: 0.37205 ~—:0.37307 
50; 0.41185 0.41135 KS, AD, AND CV SIGNIF. LEVELS NOT 
15: 0.4597 (0.45597 EXACT WITH ESTIMATED PARAMETERS. 
90: 0.49305 0.4953 
95: 0.5246 0.52209 NOTE: A SMALL SIGNIFICANCE LEVEL 


(EG. Ps.01) INDICATES LACK OF FIT 
CHI-SQUARE GOODNESS OF FIT TABLE 


LOWER UPPER OBS BP = OE (OE) 2) 
“INF. 0.28429 5.3893 3.6107 2.4191 
0.28429 0.31983 28 35.812 -7.8116 1.7039 
0.31983 0.35536 115 115.1 0.10135 0.000089245 
0.35536 0.3909 199 206.12 7.121 0.24601 
0.3909 0.42644 240 235.71 4.2893 0.078055 
0.42644 0.46197 193 190.28 2.7187 0.038843 
0.40197 0.49751 126 116.89 9.1094 0.7099 
0.49751 0.53305 57 57.85 0.85013 0.012493 
0.53305 0.56858 26 24.101 1.8989 0.14961 
0.56858 0.60412 5 8.7464 -3.7464 1.6047 
0.60412  +INF. 2 3.9966 -1.9966 0.99744 
TOTAL 1000 1000 7.9602 


«oO 


Figure 16. Lognormal Graphical Analysis of Randomly Generated Numbers 
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APPENDIX E. PROBABLISTIC TIME DISTRIBUTIONS 


The figures in this appendix depict the command and control sequence and timeline 


distributions for the digitized and non-digitized simulations. 
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Figure 17. Digitized C2 Sequence and Timeline Distribution 


73 


Transmit 
CFF from 
FO 










TARGET 
ENTERS BODE 
CELL AT S2 
NODE - 
PASSED TO 
FSE 








NO 





FSO FSO 
AVAILABLE AVAILABLE 


CUEING 
TIME 






Oo 5) 20 





YES OF Seo 
NO 
LNO ] 
CUEING 
Oo 118 20 


YES 






TRANSMIT 
FIRING ORDER 
TO BATTERY 








PLT LOR 
PREDICTION 
TIME 







PET LOR 
PROCESSING 
TIME 







PLT LOR 
ISSUES FIRING 
ORDER 





0 418 45 






DATA aad 
LAUNCHER 





Figure 18. Non-Digitized C2 Sequence and Timeline Distributions 
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APPENDIX F. INPUT DATA TEXT FILE 


Modifications to the text file shown in this appendix allow the user to change 


parameters within the model execution. 


Data Text File 


{l. Number of ave of approach into EA. USED IN ACQ TIME FUNCTION} 

= Ave height of vegetation in HRS. Used in Acquisition Time Function} 
a Maximum time in minutes FO spends SEARCHING before reinitialization} 
aa RATE OF MARCH OF ENEMY VEHICLES IN meters per second} 


S55 

{5. EFOGM VELOCITY IN meters per second} 
Bee 240 

{6. MAX EFOGM RANGE IN Kilometers} 

15.0 


{7. EXTRA TIME FOR EACH ADDITIONAL RADIO CONTACT ATTEMPT IN MINUTES} 
ee OF FO. DETECTION by -fHE BNEMY} 

eee OF FO DESTRUCTION GIVEN DETECTION BY ENEMY} 

AG SECTOR OF SCAN FOR OBSERVER/SENSORS IN DEGREES} 

a SIZE OF FIELD OF VIEW OF OBSERVER/SENSOR IN DEGREES} 

co. MEAN TIME IN MINUTES TO XMIT RADIO CALL TO BDE HQ FROM OBSERVER} 
ee SecIGMA STD DEV IN MINUTES OF TIME REOD TO XAMIT RADIO CALL TO BDE HO FROM FO} 
ee OF FIMES FO WILL ATTEMPT COMMO WITH BDE SEFORE: GIVING UF} 
fae: MEAN TIME IN MINUTES FOR BDE RTO TO PROCESS CFF FROM FO} 

alee 3 SIGMA STD DEV IN MINUTES FOR BDE RTO TO PROCESS CFF FROM FO} 
nee MEAN TIME IN MINUTES FOR RTO TO CHECK THE STATUS OF FSO} 

cr Seo lGMA STD) DEV IN MINUTES FOR RTO TO CHECK. THE STATUS OF FSO} 
0202778 

{19. MEAN TIME IN MINUTES FOR FSO TO PROCESS CFF} 


0.41667 

{20. 3 SIGMA STD DEV IN MINUTES FOR FSO TO PROCESS CFF} 
0.0611 

(ieee neeN TIME IN MINUTES TO CHECK ON EFOGM C2 CELL} 

0.25 

j2220)3 SIGMA STD DEV IN MINUTES TO CHECK ON EFOGM C2 CELL} 
0702778 


{23. MAX PROCESSING TIME IN MINUTES BEFORE REACHING THE EFOGM LNO THAT A CFF IS CANCELED 
BY FSO/LNO} 

ie. 0 

{24. TIME IN MINUTES FOR EFOGM LNO TO DETERMINE IF TARGET HAS ALREADY BEEN PROCESSED} 
oe 25 

{25. MEAN TIME IN MINUTES FOR EFOGM LNO TO PROCESS CFF} 

C2 716677 

{26. 3 SIGMA STD DEV IN MINUTES FOR EFOGM LNO TO PROCESS CFF} 

020667 

{27. MEAN TIME IN MINUTES FOR EFOGM C2 CELL TO XMIT FIRING ORDER TO BATTERY} 

0225 

{28. 3 SIGMA STD DEV IN MINUTES FOR EFOGM C2 CELL TO XMIT FIRING ORDER TO BATTERY} 


1 


0.1667 
t2e. 
U2 353 
{30 
OL 1TSG 
(Geille 
02333 
(32. 
0.0621 
(337 
G225 
{34. 
021667 
fie 
O79 66 
126: 
O20277 
{37% 
03333 
(oo. 
i 
139: 
Os 
{40. 
OZ0L 
{41° 
O05 
{42. 
O02 
{43. 
O05 
{44. 
O04 
{45. 
O02 
{46. 
0.01 
{47. 
O01 
{48. 
C265 
{49. 
0.0L 
{50. 
O02 
{fale 
6.95 
52% 
0202 
Se 
Oe OCs: 
{54. 
O02 
toe 
On 
{26- 
O.2 
tSa. 
On 
face 
0. 0s 
{39° 
0.04 
£60; 
Oo 
foi] 
O22 
toa 
0.4 
(63% 
0.02% 
{64. 
OL 


MEAN TIME IN MINUTES FOR PLT LDR TO MAKE PREDICTION} 


3 SIGMA STD DEV IN MINUTES FOR PLT LDR TO MAKE PREDICTION} 


7 


MEAN TIME IN MINUTES FOR PLT LDR TO) PROCESS CEE} 


3 SIGMA STD DEV IN MINUTES FOR PLT LDR TO PROCESS CFF) 


il 


MEAN TIME IN MINUTES FOR PLT LDR TO XMIT FIRING ORDER TO SQDS} 
3 SIGMA STD DEV FOR PLT LDR TO XMIT FIRING ORDER TO SQDS} 


MEAN TIME IN MINUTES FOR GUNNER TO APPLY DATA TO LAUNCHER} 


a 


3 SIGMA STD DEV IN MINUTES FOR GUNNER TO APPLY DATA TO LAUNCHER} 


8 


TIME IN MINUTES TO GUNNER TO EXECUTE MISFIRE PROCEDURES } 


TIME IN MINUTES FOR GUNNER TO SELECT,ARM, INITAILIZE NEW MISSILE IF DUD OCCURS} 


TIME IN MINUTES “FOR VEO 


transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 
transition 


transition 


probability 
probability 
probability 
probability 
probability 
probability 
probability 
probability 
probability 
Probability 
probability 
probability 
probability 
probability 
probability 
probability 
probability 
probability 
probability 
probability 
Probability 
probability 
probability 
probability 


probability 


TO EVALUATE BDA AND DECIDE ON COURSE OF ACTION} 


from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 
from state 


from state 


11 to state 


12 to state 


14 to state 


16 to state 


18 to state 


18 to state 


20 to state 


20 to state 


20 to state 


21 to state 


21 to state 


22 to state 


24 to state 


23 to state 


26 to state 


28 to state 


29 to state 


31 to state 


Bul state 


31 to state 


32 state 


32 state 


32 state 


30 state 


33 state 
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2 
ita, 
16} 
67 
Mee) 
20} 
19} 
20} 
21} 
1:3} 
20} 
ZS} 
Za 
ZS) 
27) 
55) 
Sy) 
36) 
5) 
34} 
36} 
esi), 
34} 
a8 


37) 


{65. transition probability from state 34 to state 37} 
{66. transition probability from state 35 to state 37} 
{67. transition probability from state 36 to state 38} 
{68. MEAN TIME IN MINUTES FOR OBSERVER TO XMIT CFF TO AFATDS} 


fee orD DEY IN MINUTES FOR OBSERVER TO XMIT CFF TO AFATDS} 
0.0278 

{70. MEAN TIME IN MINUTES FOR AFATDS TO EVAL TGT DATA} 

O.25 

(7 ieeeorD DEY IN MINUTES FOR AFATDS TO EVAL TGT DATA} 

020556 

{72. MEAN TIME IN MINUTES FOR AFATDS TO XMIT DATA TO BDE/PL/GUNS} 
Oeil 67 

(22st D DEV IN MINUTES FOR AFATDS TO XMIT DATA TO BDE/PL/GUNS} 
O201667 

{74. MEAN TIME IN MINUTES FOR BDE TO DETERMINE IF GOOD TGTING SOLN} 
Oz25 

jo otD DEV IN MINUTES FOR BDE TO DETERMINE IF GOOD TGTING SOLN} 
0.0556 

(76. MEAN TIME IN MINUTES FOR BDE TO SEND CONFIRMATION TO PL} 

O.. Tu67 

(gem otD DEV IN MINUTES FOR BDE TO SEND CONFIRMATION TO PL} 
CF01667 

{78. MEAN TIME IN MINUTES FOR PL TO SEND CONFIRMATION TO GUNS} 
C1667 

(79° STD DEV IN MINUTES FOR PL TO SEND CONFIRMATION TO GUNS} 

S. 01667 

{80. MEAN TIME IN MINUTES FOR GUNNER TO XFER DATA TO SYSTEM} 
Oe | 

hoteeolD DEV IN MINUTES FOR GUNNER TO XFER DATA TO SYSTEM} 


0.0833 

{82. MEAN TIME IN MINUTES FOR GUNNER TO SELECT AND VALIDATE ROUTE} 
i. 330 

{83. STD DEV IN MINUTES FOR GUNNER TO SELECT AND VALIDATE ROUTE} 


D2 0633 
{84. MEAN TIME IN MINUTES FOR GUNNER TO APPLY DATA TO LAUNCHER} 


O3.25 
{85. STD DEV IN MINUTES FOR GUNNER TO APPLY DATA TO LAUNCHER} 


0.01667 


7] 





APPENDIX G. TARGET DATA FILE 


The following table provides the data incorporated into the input text files. It 
provides the data only for the high value targets located within HRS 33.5 and defined by 
the DBBL. The data contain values for the ranges from each high value target to each of 
the five sensor/FOs used in the simulation, as well as the ranges to the EFOGM launcher. 
Target height is provided as an input to the DYNTACS detection algorithm (discussed in 


Appendix F). See the input text files for formatting. 


TARGET TARGET POSITION IN INITIAL INITIAL INITIAL TARGET 
NUMBER DESCRIPTION ATTACK RANGE TO RANGE TO RANGE TO HEIGHT (m) 
COLUMN SENSORS 1-3 | SENSORS 4-5 EFOGM 
km km BATTERY (km 
RECON CO 4 KM BACK IN 
ELEMENT 
RECON CO 4 KM BACK IN 
BTR60 ELEMENT 
1ST MRC 
1ST MRC 


ZSU/23-4 #4 | 600 ACK N 
saa see ee 39.60 42.60 57.60 3.75 


ede 
| 42.70 | 87.70 | 375 


SA-13/1 ON | 750m BACK IN 
CHASSIS 

SA-13/2 ON | 800 m BACK IN 
CHASSIS 













SA-13/3ON | 850m BACK IN 
ie caine 39.85 42.85 57:65 


CHASSIS 
SA-13/4 ON 900 fa 
tee yes 39.90 42.90 57.90 oe 
CHASSIS 
BCDR 
ARTY MOBILE 1250 m BACK |; 
“SMALL FRED" 
MRC 2CDR 1800 m BACK 
16 BG BS omIEAC | 4070 43.70 | 58.70 | 70 
MRC 3 CD 
a eC Ame ates | 4435 oe 35 
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ok per’ ook ae 
ee) NO — 








TARGET TARGET POSITION IN INITIAL INITIAL INITIAL TARGET 
NUMBER DESCRIPTION ATTACK RANGE TO RANGE TO RANGE TO HEIGHT (m) 
COLUMN SENSORS 1-3 | SENSORS 4-5 EFOGM 
(km) km BATTERY (km 
1st D-30 HOW | 2900 m BACK 
BTRY FDC - IN1ISTMRB | 41.80 44.80 99.80 
BTR 60 
D-30 HOW BN | 3250 m BACK 
Toocerheo| Nistems | 2215. | oor 
2nd D-30 HOW | 3300 m BACK 
pre a Rapa eae ! 42.20 45.20 60.20 2.31 
BTR 60 


3rd D-30 HOW | 3650 m BACK 
21 UD SO HOW, | SSs0m BACK |g? Souler 60.55 2.31 
BTR 60 


fy pins a atl teas aco | ghia | 


= eee er ee Pee [ar 
"SMALL FRED" 
or eee eto 


D-30 HOW BN } 2750m BACK 
| TOC - BTR 60 IN 2ND MRB 


BTRY FDC - 
BTR 60 


3rd D-30 HOW | 3150m BACK 
BTR 60 
ZPU-4 /#1 FRONT OF 
GROUP 
ZPU-4/#2 50 m IN BACK 
GROUP 
ZPU-4/#3 100 m IN BACK 
GRP 
ZPU-4/#4 150 m IN BACK 
GRP 
ZPU-4/#5 200 m IN BACK 
GRP 
ZPU-4/#6 250 m IN BACK 
: 36 | pO MIN BAC 50.05 5@.05 68.05 a 
— rong : 
277 


M1939 300 m IN BACK | 4 | 
37mm(SP) ADA | OF BDE HQ 90.10 53.10 | 68.10 


GUN #1 GRP 


M1939 350 m IN BACK 
37mm(SP) ADA OF BDE HQ 90.15 93.15 68.15 
GUN #2 GRP 
M1939 400 m IN BACK 
GUN #3 GRE 
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TARGET TARGET POSITION IN | INITIAL INITIAL INITIAL TARGET 
NUMBER DESCRIPTION ATTACK | RANGE TO RANGE TO RANGE TO HEIGHT (m) 
COLUMN | SENSORS 1-3 | SENSORS 4-5 EFOGM 
| km) km BATTERY (km 
ey 450 m IN BACK 
pe #4 GRP 
|: CDR BMP } 500 m IN BACK 
GRP 
BDE S38 BMP | 550m IN BACK 
GRP 
ARTY MOBILE | 800 m IN BACK 
AT MOREE [PERNA] 50.60 | 5860 | 68.60 
"SMALL FRED" GRP 
ZU-23 (T) ADA | 950 m IN BACK 
GRP 
ZU-23 (T) ADA 1000 m IN 
HQ GRP 
ZU-23 (T) ADA 1050 m IN 
HQ GRP 


ZU-23 (T) ADA OOmIN | 
any ie oreee | 20:98 93.90 68.90 


HQ GRP 


1st M46 130mm 2100 m IN 
FDC BTR 60 HQ GRP 
M46 130mm 2450 m IN 
TOC BTR 60 HO GRE 
2nd M46 2500 m IN 
130mm BACK OF BDE 92.30 99.30 70.30 
HOW(T) BTRY HQ GRP 
FDC BTR 60 

3rd M46 130mm 2850 m IN 
FDC BTR 60 HQ GRP 
4th ae 130mm 3250 m IN 
oe BTR 60 HQ GRP 
M46 130mm 3600 m IN 

HOWCT) 2nd | BACK OF BDE 53.40 56.40 71.40 Zon 
BN TOC BTR HQ GRP 

60 

5th M46 130mm 3650 m IN 
FDC BTR 60 HQ GRP 
6th M46 130mm 4000 m IN 
FDC BTR 60 HQ GRP 
ist BM-21 4400 m IN 
HQ GRP 





BM-21BNTOC| 5750 mIN 
By arate aren aes Sy he s) t2-50 


HQ GRP 
2ND BM-21 5800 m IN 
BTRY FDC BACK oF BDE 94.60 97.60 72.60 
HQ GRP 
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TARGET TARGET POSITION IN INITIAL INITIAL INITIAL TARGET 
| NUMBER DESCRIPTION ATTACK RANGE TO RANGE TO RANGE TO HEIGHT (m) 
| COLUMN SENSORS 1-3 | SENSORS 4-5 EFOGM 
(km) | km BATTERY (km 
3RD BM-21 6150 m IN 
as Oe ae 54.90 | 57.90 ee 2.89 
HQ GRP 
a CO CDR #1 | 150 m IN BACK 
eet QMINBACK) 57.45 | || WeOoueeey odo Pe a 
FRONT 
TK BN see T- | 650m IN BACK 
FRONT 
ARTY MOBILE | 700m IN BACK | 7, 
RECON POST OF TK BN 08.00 61.00 6.00 
"SMALL FRED" FRONT 



































TK CO CDR #2 1000m IN 
ae ae SOON a oere5 61.25 76.25 
BN FRONT 
TK CO CDR #3 1650 m IN 
BN FRONT 
200 m IN BACK 
OF MRB 79. 00 
* FRONT 
66 | MRB CDR BMP | 650 m IN BACK 
OF MRB 
FRONT 





ARTY MOBILE | 750 
"SMALL FRED" FRONT 
MRC 2 CDR 1300 m IN 
FRONT 
MRC 3 CDR 1950 m IN 
FRONT 
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APPENDIX H. THE DYNTACS MODEL’ 


Any discussion of the DYNTACS equation must first consider the history of 
detection modeling. B.O. Koopman first published a report on target acquisition modeling 
for the U.S. Navy in 1946. In his report he defined detection as “...that event constituted 
by the observer’s becoming aware of the presence and possibility of the position and even 
in some cases the motion of the target.” [Ref. 2, pg. 67] There are several distinct areas 
of target acquisition. These include: 


e Cueing Information, which provides the approximate area on the battlefield for 
continuing search (1.e., an explosion, gun flash, reflections). 


e Detection means that an observer has decided that something within his field of 
view has some military value (e.g. he determines it’s a man and not a deer). 


e Classification occurs when the observer is able to distinguish between the type of 
target discovered (e.g. a wheeled vehicle versus a tracked vehicle). 


e Recognition occurs when the observer 1s able to determine between a smaller 
variety of tracked vehicles (e.g. a tracked troop carrier versus a tank). 


e Identification occurs when the observer is able to precisely identify the target (e.g., 
a BRDM versus a BTR-60P). 


The continuous looking model was the second of Koopman’s basic detection 
model paradigms. The continuous looking model is based on a detection rate function 
D(t), which is the probability of detecting a target in the short time interval t ,which is 
proportional to the length, A, of the time interval. The formula for this detection function 


is: 


P(detect in [t, t + AT]) = D(t) *AT. (1) 


' This discussion of the DYNTACS equation is paraphrased from Ref. 2, Chapter IV. 
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Assuming that the Detection function remains constant for all values of t, then for a longer 
period of time T, T = N * AT, where N is the number of intermittent opportunities to 
detect a target, or “glimpses”. Therefore, the probability of detection in a time period of 
length T is: 
P(Detect in length T) = I - P(fail to detect in N tries) 

= {| -(1 -D*AT)” 

= (/8@=D*T/N) (2) 
Taking the limit, as N approaches infinity, AT approaches 0, and T = N*AT held constant, 


this becomes 


P(Detect in length T) = I - EXP(-D*T). ) 


This equation is recognized as the cumulative distribution function (CDF) of the 
exponential probability distribution. This detection rate function can be used in both fixed 
time and event scheduled simulations. 

The DYNTACS curve fit model was developed in conjunction with a set of field 
experiments conducted in the 1960’s to provide good experimental data for detection 
rates as a function of several parameters. A “detection” within these experiments 1s the 
same as an “identification” listed above. The experiment has some limitations. 


e Allexperiments were conducted during daylight at limited ranges (1.5 
kilometers). 


e Observers had restricted fields of view (30 degrees). 

e Observers were stationary, and not mounted in vehicles. 

e Observers used the naked eye for observation with no vision aids. 

The experiments verified that the exponential distribution could be used for 
observed detection times, stationary targets, and that a constant detection rates were valid. 


The experimental data were used in regression models to compute detection rates D for 


various situations. 
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While the resulting detection rate functions is only valid for simulations which 
follow the experimental situation, it has been extended to ranges of 5 kilometers and for 


observers from vehicles. Using Equation (3), and solving for t yields: 
t = In (1- P(Detection))/-D. (4) 


The value for P(Detection) can easily be simulated using a Monte Carlo draw from a 
generator of Uniform(0,1) random variables. The DYNTACS experiment provides the 


value of D, where 
D = P,, * {-0.003 + [ 1.088 / DENOMINATOR] }, (>) 
and 
DENOMINATOR = 1.453 + *(0.05978 + 2.188R° - 0.5038*CV). (6) 
Observation conditions for DENOMINATOR are given by 
e t= terrain complexity code. This is equivalent to the number of potential 


avenues of approach for the enemy into the battle area. 


e R= apparent range in kilometers. The apparent range is the range at which the 
image of a fully exposed M60 tank would be the same height as an image for 
the current target. R = (actual range * M60 height) / (target height * percent 
visible). Percent visible is the quantity target height minus the average 
vegetation height in the area, divided by the target height. 


e CV -=crossing velocity in meters per second. This is the perpendicular 
component of the target velocity relative to the observer target line. 


e P, =the probability the observer is looking in the 30 degree search sector 
which contains the target. 


Target acquisition within the simulation model by the function TimeToAcquire 
which uses the DYNTACS equation. TimeToAcquire is called by the sub-routine Search 
within the simulation model sub-programs MASTERTIMER and MASTERTIMER2 The 
values above are provided by the user within the input data file. Once a value for D is 


calculated by TimeToAcquire, this value is substituted into equation (4) within the 


function, which calculates the time of detection for that specific target. Once Search has 
looped over all possible targets for a specific FO, the detected target with the minimum 
time of acquisition becomes the acquired target for that FO, and is forwarded for further 


processing by the simulation model. 
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APPENDIX I. SAMPLE OUTPUT DATA FILE 


Table 6 below contains a portion of an output data file for a non-digitized 


replication, after formatting in a spreadsheet program. 














LAL/ 
DUD 
Missile 
CFF CFF Time 

|(O,T)| State 1 | State2 | State3| State 4-| State5 | State6 | State7 | State8 | State 9 | State 10 | State 11 | 
(1, 1)| 24.695 | 24.985 | 26.264 | 9099 | 9999 | 99099 | 9999 | 9999 | 9999 | 9999 | 990.9 | 
1,2 
(1, 3)| 92.422 | 92.655 | 93.693 | 93.846 | 94.574 | 94.795 | 96.18 | 98.693 | 999.9 | 122.165 | 124.42 _ 
(1, 4)| 94.782 | 95.052 | 96.013 | 96.013 | 96.638 | 97.276 | 98.111 | 100.604 | 999.9 | 122.668 | 124.90_ 
(1,5)| 97.76 | 98.018 | 98.949 | 99.992 | 100.574 | 102.925 | 104.3 | 106.551 | 999.9 | 122.726 | 124.97 | 
if1,/6) | 109.912 | 110.124 | 111.101] 9099 | 9999 | 9999 | 9999 | 9999 | 999.9 | 9999 | 999.9 | 
(1, 9) | 107.302 | 107.549 | 108.457| 999.9 | 999.9 | 999.9 | 999.9 | 999.9 | 999.9 | 999.9 | 999.9 | 
(1,10) 90.042 | 90.267 | 91.314 | 91.314 | 92.046 | 92.046 | 93.187 | 95.338_| 999.9 | 124.511 | 126.73 | 
(1,11)| 106.211 | 106.44 | 107.295} 107.367 | 107.899 | 110.933 | 111.982 | 114.235 | 999.9 | 124.226 | 126.47 | 
|(1,12)| 118.809 | 119.02 | 120.202] 999.9 | 999.9 | 999.9 | 999.9 | 9999 | 9999 | 999.9 | 999.9 


Table 6. Example Output Data File from Non-Digitized Replication Run 


CFF | LNO 
arrives @| starts |arrives to] Processing 
LNO [Processing Battery | Complete 


Detection CFF 
Time _ |arrives @| arrives 
BDE @ FSO | processing 



































Output tables similar to Table 6 were produced for each replication of the 
simulation model. From these tables, the time a CFF from a specific FO on a specific 
target (enumerated in each row) arrive al various processing stages can be seen. “999.9” 
is the default value that the data array OutData (discussed in Chapter 3) is initialized to at 
the beginning of each replication. For example, the CFF from FO 1 on Target 2 ( shown 
in row four) indicates that Target 2 was detected by FO 1 (seen under State 1) at 12.5 
minutes into the battle. The CFF reached BDE (time under State 2) at 12.69 minutes in 
the battle, and that the FSO received the CFF (time under State 3) at 13.72 minutes into 
the battle. From that point, all remaining column values are “999.9”, indicating that the 
CFF never went beyond the FSO. This could have been because the total CFF processing 
time was too great, or that the FSO had previously processed a CFF on that target from 


another FO. 
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In order to obtain total CFF processing time, values in State 8 (Battery completes 
CFF processing) are subtracted from values in State 1 (Detection Time),and the total CFF 
processing time 1s obtained. Likewise, by subtracting the time under State 4 from State 3, 
or the time under State 6 from State 5, the time a CFF spent in the FSO or LNO queue 
can be found. Similar processing of all CFF was conducted for both digitized and non- 
digitized replications. These data values were then used to compute the statistics for CFF 
processing time and LNO and FSO queue time, which are used in this thesis. A summary 


of all CFF processing times is located in Appendix N. 
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APPENDIX J. LOGNORMAL DISTRIBUTION FIT, NON-DIGITIZED, BASE 
CASE DATA 


The following figures depict the AGSS analysis of the non-digitized, base case data 


to the lognormal distribution. 


ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA : DIGIT 

SELECTION : ALL 

X AXIS LABEL: DIGIT 

SAMPLE SIZE : 1200 

CENSORING : NONE 

FREQUENCIES : 1 

EST. METHOO : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 


CONF. INTERVALS = COVARIANCE MATRIX OF 
E (95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER MJ S| GMA 
MJ 0.37196 0.36692 0.37699 6 .S905E-6 0.0000E0 
S|GMA 0.088931 0.085545 0.092678 0.0000EO 3.2953E-4 
LOG LIKELIHOOD FUNCTION AT MLE = 754.81 


SAMPLE FITTED GOODNESS OF FIT TESTS 
MEAN =: 1.4564 1.4563 CHI-SQUARE : 26.665 
STO DEV : 0.1326 0.12977 DEG FREED: 6 
SKEWNESS: 0.82888 0.26803 SIGNIF : 0.00016727 
KURTOSIS: 5.1548 3.128 KOLMESMIRN : 0.028479 
+ BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF : 0.28469 
CRAMER-V M: 0.16116 
PERCENTILES SAMPLE —-FITTED SIGNIF =: > .15 
5: 1.2615 253° ANDER-DARL : 1.2424 
1.3015 2943 SIGNIF =: > 15 
1.364 3662 
1.447 4506 KS, AD, AND CV SIGNIF. LEVELS NOT 
1.5345 5402 EXACT WITH ESTIMATED PARAMETERS. 
1.619 6257 
1.6795 6791 NOTE: A SMALL SIGNIFICANCE LEVEL 
(EG. Ps.01) INDICATES LACK OF FIT 


CHI-SQUARE GOOONESS OF FIT TABLE 


UPPER 08S EXP O-£  ((O-£)*2)=€ 
1.1895 5 15.40! -10.401 7.0241 
281 77 81.89 4.8899 0.29199 
13725 247 223.05 23.952 --2.572 
464 333 329.19 3.806 0.044004 
5555 295 291.12 3.8782 0.051665 
167.39 11.386 0.77452 
66.913 -10.913 1.7798 
19.659 -2.6588 0.3596 
+INF. 14 5.3876 8.6125 13.768 
1200 1200 26.665 





Figure 19. Lognormal Distribution Graphical Analysis 
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DIGITIZED DATA EDA - LOGNORMAL DISTRIBUTION 


LOGNORMAL DENSITY FUNCTION, N=1200 LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=1200 
S 
ES et 


0.20 0.25 
0.6 0.8 


RELATIVE FREQUENCY 
0.10 0.16 


CUMULATIVE PROBABILITY 
0.4 


0.068 
0.2 





i; f 
al j 


1.2 14 1.6 1.6 2.0 2.0 
DIGIT DIGIT 


LOGNORMAL PROBABILITY PLOT, N=1200 LOGNORMAL CUMULATIVE HAZARD FUNCTION, N=1200 





12 14 1.6 16 2.0 1.2 1.4 1.6 
DIGIT DIGIT 


Figure 20. Lognormal Distribution Data Analysis 
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APPENDIX K. WEIBULL AND BETA DISTRIBUTION COMPARISONS, NON- 
DIGITIZED, BASE CASE + 25% DATA 


The following figures depict the AGSS analysis of the non-digitized, base case + 


25% data for the Weibull and Beta distributions. 


ANALYSIS OF WEIBULL OISTRIBUTION FIT 


OATA 

SELECTION : 

X AXIS LABEL: 

SAMPLE SIZE : 

CENSORING =: 

FREQUENCIES : 

EST. METHOD : MAXIMUM LIKELIHOOD 

CONF METHOD : ASYMPTOTIC NORMAL APPROXIMATION 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER c a 
C (SHAPE) 2.7083 2.5859 2.8307  0.0038977 0.0042827 
a (SCALE) 19.292 18.867 19.717 0.0042827 0.047025 
LOG LIKELIHOOD FUNCTION AT MLE = -3993.4 


SAMPLE FITTED GOOONESS OF FIT TESTS 
MEAN ame (Fie (OF, 17,158 CHI-SQUARE : 61.689 
STD OEV: 6.9463 6.8343 OEG FREED: 10 
SKEWNESS: 0.16079 0.64783 SIGNIF : 1.7336E-9 
KURTOSIS: 2.0611 1.2917 KOLM-SMIRN : 0.053485 
s BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF =: 0.0020983 
CRAMER-V M: 0.84601 
PERCENTILES SAMPLE FITTED SIGNIF :  < .01 
S: : ANDER-DARL : 6.8602 
SIGNIF : ¢.01 


KS, AD, AND CV SIGNIF. LEVELS NOT 
EXACT WITH ESTIMATED PARAMETERS. 


NOTE: A SMALL SIGNIFICANCE LEVEL 
(EG. Ps.01) INDICATES LACK OF FIT 


CHI-SQUARE GOODNESS OF FIT TABLE 


UPPER OBS EXP O£  ((O-£)e2)2 
7.4723 127 88.448 38.552 16.804 
9.963 101 95.981 5.0193 0.26248 
12.454 124 131.34 -7.3403 0.41023 
14.945 130 156.6 -26.595 4.5166 
17.435 152 166.08 -14.078 1.1933 
19.926 130 158.07 28.065 4.9831 
22.417 115 135.42 -20.417 3.0782 
24.908 136 104.42 31.579 9.5503 
27.398 92 /2.338 19.662 5.3442 
29.889 62 44.888 17.112 6.5235 
32.38 22 24.859 -2.8594 0.32889 
34.871 7 12.237 5.2372 2.2414 
INF. 1 B.3334 -7.3334 6.4534 
1199 1199 61.689 





Figure 21. Weibull Distribution Data Analysis 
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NON-DIGITIZED DATA EXCURSION ONE/ WEIBULL EDA 


WEIBULL DENSITY FUNCTION, N=1199 


0.08 


RELATIVE FREQUENCY 


0.04 


0 10 20 30 
NONDIGIT 


WEIBULL PROBABILITY PLOT, N=1199 





2 A 10 20 40 
NONDIGIT 


WEIBULL CUMULATIVE DISTRIBUTION FUNCTION, N=1199 





2 
3 
2 
Be 
g 6 
py 
x 
Eo 
: 
. 

6 

: 

10 20 30 
NONDIGIT 


WEIBULL CUMULATIVE HAZARD FUNCTION, N=1199 


6 


«b 


CUMULATIVE HAZARD 


\ 


anes 
10 20 30 
NONDIGIT 


Figure 22. Weibull Distribution Graphical Analysis 
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ANALYSIS OF BETA DISTRIBUTION FIT 


DATA > NORMNONDIGIT 

SELECTION : ALL 

X AXIS LABEL: NORMNONDIGIT 

SAMPLE SIZE : 1199 

CENSORING =: NONE 

FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 

CONF WETHOD : ASYMPTOTIC NORMAL APPROXIMATION 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) _PARAVETER ESTIMATES 
PARAMETER ESTIMATE LOWER UPPER Q 


P 2.1297 2.5199 2.9396 0.011457 0.010378 
Q 2.9243 2.6982 3.1504 0.010578 0.013307 
LOG LIKELIHOOD FUNCTION AT MLE = 295.01 
SAMPLE —- FITTED GOODNESS OF FIT TESTS 
MEAN =: (0.47874 0.4828 CHI-SQUARE : 19.458 
STD DEV: 0.19439 0.19372 DEG FREED: 9 
SKEWNESS: 0.16079 0.046415 SIGNIF =: 0.021566 
KURTOSIS: 2.0611 2.3095 KOLM-SMIRN : 0.046718 
+ BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF : 0.010666 
CRAMER-V M : 0.4386 
PERCENTILES SAMPLE = FITTED SIGNIF : < .10 
9: 0.18328 §=— 0. 16823 ANDER-DARL : 4.1562 
10: 0.20446 = 0.22422 SIGNIF =: < .0f 
2: 0.30929 0.33691 
90: 0.47225 = 0.48067 KS, AD, AND CV SIGNIF. LEVELS NOT 
3: 0.65576 0.62662 EXACT WITH ESTIMATED PARAMETERS. 
90: 0.75076 0.74466 
99: 0.79929 0.80478 NOTE: A SMALL SIGNIFICANCE LEVEL 


(EG. Ps.01) INDICATES LACK OF FIT 
CHI-SQUARE GOODNESS OF FIT TABLE 


LOHER UPPER 08S EXP OE (OE) v2) 38 
-INF. 0.20911 127 101.6 25.307 6.3483 
0.20911 0.27862 101 96.561 4.4395 0.20411 
0.27882 0.34852 124 123.75 0.25514 0.00052607 
0.34852 0.41822 130 142.46 12.457 1.0893 
0.41822 0.48793 152 150.84 1.1559 0.0086572 
0.48793 0.55763 130 149.06 -19.063 2.4378 
0.55763 0.62733 115 136.29 -21.293 3.3265 
0.62733 0.69704 136 115.26 20.739 3.7316 
0.69704 0.76674 92 87.934 4.066 0.18801 
0.76674 0.83644 62 57.803 4.197 0.30474 
0.83644 0.90615 22 29.294 -7.2936 1.816 
0.90615 +INF. 8 ~—«8.1435. 0.14352 0.025292 
TOTAL 1199 1199 19.458 


Figure 23. Beta Distribution Data Analysis 
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NON-DIGITIZED DATA EXCURSION ONE/BETA EDA 





BETA DENSITY FUNCTION, N=1199 BETA CUMULATIVE DISTRIBUTION FUNCTION, N=1199 
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Figure 24. Beta Distribution, Graphical Analysis 


94 


APPENDIX L. LOGNORMAL DISTRIBUTION COMPARISON, DIGITIZED, 
BASE CASE DATA 


The following figures depict the AGSS analysis of the digitized, base case data for 


the lognormal distribution. 


ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DIGIT 
ALL 
DIGIT 


DATA 
SELECTION 

X AXIS LABEL: 
SAMPLE SIZE : 1200 
CENSORING : NONE 

FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 


COVARIANCE MATRIX OF 

PARAMETER ESTIMATES 
MU S|GMA 
6.5905E-6 0.0000E0 
0.000060 3.2953E-6 


CONF. INTERVALS 
(95 PERCENT) 

PARAMETER ESTIMATE LOWER UPPER 

MW 0.37196 0.36692 0.37699 

SIGMA —- 0.088931 0.085545 0.092678 

LOC LIKELIHOOD FUNCTION AT MLE = 754.81 


MEAN 


STD DEV : 
SKEWNESS: 0.82888 
KURTOSIS: 5.1548 


+ BASED ON MIDPOINTS OF FINITE INTERVALS 


PERCENTILES 


5: 


UPPER OBS = EXP 
sWGd0~ 0 
. 281 
3725 247 
.464 
f0de2 4290 
647 
57569 
.83 
+INF. 


eee eee 


SAMPLE 
1.4564 
0.1326 


7] 
533 


156 
96 
17 
14 


1200 1200 


SAMPLE 


FITTED GOODNESS OF FIT TESTS 
1.4563 > 26.665 
0.12977 ee 
0.26803 0.00016727 
3.128 0.028479 

0.28469 

0.16116 
Polo 
1.2424 
>=..10 


SIGNIF : 
ANDER-DARL : 
SIGNIF 


FITTED 
2008 
o20t0 
. 3662 
4506 
9402 
6257 
.6791 


KS, AD, AND CV SIGNIF. LEVELS NOT 
EXACT WITH ESTIMATED PARAMETERS. 


NOTE: A SMALL SIGNIFICANCE LEVEL 
(EG. Ps.01) INDICATES LACK OF FIT 


CHI-SQUARE GOODNESS OF FIT TABLE 


O-£ 

-10.401 

~4 8899 
25,952 
3.806 


((0-£)+2)=£ 
7.0241 
0.29199 
2.572 
0.044004 
0.051665 
0.77452 
1.7798 
0.3596 

13.768 

26. 665 


15.401 
81.89 
223.05 
529.19 
291.12 3.8782 
167-39). =) 000 

66.915 ~-10.913 

19.659 2.6588 

9.5876 8.6125 





Figure 25. Lognormal Distribution Data Analysis 
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DIGITIZED DATA EDA - LOGNORMAL DISTRIBUTION 
LOGNORMAL DENSITY FUNCTION, WV=1200 LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=1200 
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Figure 26. Lognormal Distribution Graphical Analysis 
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APPENDIX M. LOGNORMAL DISTRIBUTION COMPARISON, DIGITIZED, 
BASE CASE + 25% DATA 


The following figures depict the AGSS analysis of the digitized, base case + 25% 


data for the lognormal distribution. 


ANALYSIS OF LOGNORMAL DISTRIBUTION FIT 


DATA : DIGIT 

SELECTION : ALL 

X AXIS LABEL: DIGIT 

SAMPLE SIZE : 1200 

CENSORING =: NONE 

FREQUENCIES : 1 

EST. METHOD : MAXIMUM LIKELIHOOD 
CONF METHOD : EXACT 


CONF. INTERVALS COVARIANCE MATRIX OF 
(95 PERCENT) PARAMETER ESTIMATES 
PARAMETER ESTIMATE LOWER = UPPER MU SIGMA 
MJ 0.5554 0.55096 0.55984 9. 1182E-6 0.0000E0 
SIGMA 0.07837 0.075586 0.081672 0.0000E0 2.5591E-6 
LOG LIKELIHOOD FUNCTION AT MLE = 686.38 


SAMPLE FITTED GOODNESS OF FIT TESTS 
MEAN =: 1.748 1.748 CHI-SQUARE : 2.3033 
STO DEV: 0.13836 Oi 1372 DEG FREED: 
SKEWNESS: 0.50372 0.23596 SIGNIF : 0.80579 
KURTOSIS: 4.6185 5.0991 KOLM-SMIRN : 0.010715 
« BASED ON MIDPOINTS OF FINITE INTERVALS SIGNIF <: 0.99915 
CRAMER-V M: 0.016625 
PERCENTILES SAMPLE FITTED SIGNIF == > 19 
2: 1.005 Losts ANDER-DARL : 0.17351 
125/59 1.5761 SIGNIF =: > .15 
1.6535 1.653 
1.7425 1.7426 KS, AD, AND CV SIGNIF. LEVELS NOT 
1.835 1.8372 EXACT WITH ESTIMATED PARAMETERS. 
1.92 1.9268 
1.975 1.9825 NOTE: A SMALL SIGNIFICANCE LEVEL 
(EG. Ps.01) INDICATES LACK OF FIT 


CHI-SQUARE GOODNESS OF FIT TABLE 


LOWER UPPER EXP OF  ((O-£)s2)+£ 
-INF. 1.498 32.151 4.151 0.53594 
1.498 1.605 144.13 5.8722 0.23925 
1.605 1.712 316.3 8.2956 0.21757 
1.712 1.819 356.9 8.0957 0.18363 
1.819 1.926 229.48 5.525 0.13302 
1.926 2.033 91.502 -3.5016 0.134 
2.033 2.14 24.284 4.2844 0.75587 
2.14 +IKF- 5.2603 0.73967 0.10401 
TOTAL 1200 1200 2.3033 





Figure 27. Lognormal Distribution Data Analysis 
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DIGITIZED DATA EXCURSION ONE/LOGNORMAL EDA 


LOGNORMAL DENSITY FUNCTION, N=1200 LOGNORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=1200 
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Figure 28. Lognormal Distribution Graphical Analysis 
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APPENDIX N. RUNS MATRIX WITH RESULTS 


ej 
© 


* 


in 
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